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We study the phase diagrams of a family of 3D "Walker- Wang" type lattice models, which are not 
topologically ordered but have deconfined anyonic excitations confined to their surfaces. We add a 
perturbation (analogous to that which drives the confining transition in Zp lattice gauge theories) to 
the Walker- Wang Hamiltonians, driving a transition in which all or some of the variables associated 
with the loop gas or string-net ground states of these models become confined. We show that in 
many cases the location and nature of the phase transitions involved is exactly that of a generalized 
Zp lattice gauge theory, and use this to deduce the basic structure of the phase diagram. We further 
show that the relationship between the phases on opposite sides of the transition is fundamentally 
different than in conventional gauge theories: in the Walker- Wang case, the number of species of 
excitations that are deconfined in the bulk can increase across a transition that confines only some 
of the species of loops or string-nets. The analogue of the confining transition in the Walker- Wang 
models can therefore lead to bulk deconfinement and topological order. 



I. INTRODUCTION 

Critical phenomena beyond the standard Landau- 
Ginzburg paradigm have been an active area of interest 
in condensed matter physics over the past decade. One 
area that has drawn significant interest of late is the sub- 
ject of transitions between phases of different topological 
order.f™51 

By definition, topological order cannot be diagnosed 
by any local order parameter'^. Typical hallmarks of 
this type of order include a finite ground-state degen- 
eracy that changes with the topology of the mani- 
fold, topological entanglement entropjW, and excitations 
with anyonic statistics. For phase transitions that al- 
ter the topological order, there is no known analogue of 
the Landau-Ginzburg agproach that correctly describes 
the critical properties.l^ However, in 2D there are a 
number of example s where these transitions are well 
understoocPmHUllHIll, Here we will extend this body 
of work to investigate the phase diagram of a family of 
3D modelJl^ with interesting topological properties. 

In discussing the phase diagram of a topologically or- 
dered system, there are two elements of interest: the na- 
ture of the phase transitions, and the topological orders 
of the phases themselves. Early progress on the former 
front was made in Refs. d6j and il7i where, amongst other 
things, it was shown that the confining transition of Z2 
gauge theory is in the Ising universality class. (A body of 
later numerical work has followed up on this resullPlSllSl^ 
confirming their analytical description of the phase dia- 
gram). Much later, Refs. |3]and^used similar methods to 
show that a large family of transitions in 2D Levin- Wen 
model^2SI fall into the Ising or Potts universality clas j^. 

On the second front, Slingerland and Baia^^Hll pro- 
vided a comprehensive criterion - which we use through- 
out this work - for which topological orders can be con- 
nected by a (potentially second-order) phase transition: 
two phases can be "naturally" connected by a direct 
phase transition if one of the phases can arise as the re- 



sult of forming a condensate in the other. Two familiar 
examples are the Higgs transition (Bose condensation of 
charges), and the confining transitions (Bose condensa- 
tion of vortices or vortex loops) familiar in gauge theo- 
ries. Beyond this, however, their method gives a useful 
paradigm in 2D for studying condensation transitions in 
any anyon model. It is worth emphasizing that the prob- 
lem is not trivial: the low-energy excitations in the con- 
densed phase need not be simply a subset of those of the 
uncondensed phase. 

In the present work, we will address both the question 
of the nature of the phase transition, and the relationship 
between the topological orders before and after conden- 
sation, in a family of 3D Walker- Wan^i^ type models. 
These models are not topologically ordered in the con- 
ventional sens^^H on a system with no boundaries, they 
always have a unique ground state, and admit no de- 
confined point-like excitations. However, on a 3D system 
with a 2D boundary, they exhibit surface states with chi- 
ral topological order (analogous to that of a fractional 
quantum Hall system) and deconfined anyonic excita- 
tions. Commensurate with these topological character- 
istics, their ground state wave-functions are 3D cousins 
of the loop gas or string-net ground states of Levin- Wen 
models. We refer to this set of characteristics as sur- 
face topological order, to contrast it with bulk topological 
order, in which there is a topological ground-state degen- 
eracy on a system with periodic boundary conditions (i.e. 
in the absence of surfaces). This distinction is important: 
unlike bulk topological order, a system's surface topolog- 
ical order can change without a bulk phase transition. 

Though the methodology we use here is similar to that 
employed in Ref. 3 to study condensation transitions in 
Levin- Wen models (the 2D cousins of our 3D models), 
there are some striking differences between the 2D and 
3D systems. Notably, in many cases these models admit 
condensation transitions from phases with surface topo- 
logical order to phases with 3D bulk topological order. 
This is in stark contrast to the situation in 3D gauge 



theories: the transitions we describe are analogous to 
confining transitions, involving the proliferation of vor- 
tex loops; in a gauge theory these can only decrease the 
topological ground-state degeneracy, and can never pro- 
duce a phase with bulk topological order from one with 
none. (This is because the Higgs transition necessarily 
reduces the size of the gauge group; for a gapped system 
the topological ground-state degeneracy vanishes only if 
the gauge group is trivial). 

Our key results are as follows. We identify a fam- 
ily of condensation transitions which have the same 
low energy theory as confining transitions in a Potts 
gauge theo ry. T hese transitions have b een studied both 
nume 

ricalljl23124| and theoreticallji^IMl, and are generi- 
cally first-order. Our main interest is in exposing the re- 
lationship between the pairs of 3D phases related by con- 
densation. We give examples of condensation transitions 
between Walker- Wang models of trivial bulk topological 
order and 1) the trivial phase, which has neither bulk nor 
surface topological order; 2) another Walker- Wang model 
with only surface topological order; 3) a discrete gauge 
theory, with only bulk topological order; and 4) a phase 
with both bulk and surface topological order. We will 
show that whenever bulk topological order emerges, the 
phase also has deconfined point-like excitations, which 
can be either bosons or fermions, depending on the na- 
ture of the condensed phase. In contrast, in scenarios 
2) and 4), we find that if the lattice has boundaries, the 
models have chiral anyonic excitations confined to these 
surfaces. Scenario 4) thus describes a phase with bosonic 
or fermionic point particles deconfined in the bulk, as 
well as anyons confined to the surface. 

We begin in Sect. [ll] by introducing the Potts gauge 
theory, which is a useful stepping stone for understanding 
the Walker- Wang models which we introduce in Sect. |III| 
We first discuss condensation transitions in the Abelian 
Walker- Wang models, whose Hamiltonians differ from 
those of the Potts gauge models only by additional phase 
factors. These models admit transitions in all of the cases 
(1-4) above with a minimum of complexity. 

Having understood the phase diagram of these Abelian 
models, we turn our attention to the non-Abelian case. 
We begin in Sect. [V] by discussing a family of mod- 
els with the surface topological order of SU(2)fe Chern- 
Simons theories. We show that these exhibit transitions 
analogous to those of an Ising gauge theory, as well as 
other types of transitions which do not have such simple 



analogues. In Sect. VI we describe in more technical lan- 
guage the status of transitions in general Walker- Wang 
models, giving a general prescription for determining the 
topological order of the condensed phase. 



II. REVIEW: CONDENSATION IN ABELIAN 
LATTICE GAUGE THEORIES 

We begin by examining condensation transitions in 
'Potts gauge theories I^S! ( Zp"**** models, for short), which 
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FIG. 1. Action of the plaquette operator on a plaquette 
in tlie xy-plane. On edges 1 and 2, where the arrow points 
along the x and y directions respectively, we act with We,+; 
on edges 3 and 4, where the arrow points along —x —y, we 
act with We -. 



are variations of the well-known discrete Zp gauge the- 
ories, and which give a natural point of departure from 
which to discuss transitions in the related Walker- Wang 
models.'^- We first review the models and their essential 
features, and then discuss a simple condensation tran- 
sition which will be our point of reference for under- 
standing all condensation transitions within these mod- 
els. We then examine how the topological order changes 
after the condensation. Expert readers may wish to skim 
this section for our notation; we will rely heavily on the 
framework outlined in this section both to understand 
the Walker- Wang Hamiltonians in later sections, and to 
describe the transitions in the resulting phase diagrams. 



A. The Zp 'Potts gauge theory' Hamiltonians 

The models we describe are a generalization of the 3D 
Ising gauge theorjIlSl, or 3D Toric Cod^^, to Potts spins. 
The Hilbert space of our model thus consists of a p-state 
system (a Potts spin s) on every edge of a 3D cubic lat- 
tice. To define the Hamiltonian we require two types of 
operators which act on each edge: S, which measures the 
spin, and W, which raises it. We will generally work in 
the spin basis, where a state in the Hilbert space is spec- 
ified by the value of the spin on each edge. In this basis, 
on each edge we require the four unitary operators 



Se,±|Se 

We^± |Se) = |Se ± 1 ( mod p)) 

which obey 



0,l,...p-l 



(1) 



Sei, + We2^+ = Exp [i2Tr/p Sei,e2] ^e2.+ ^ei,+ 
Sei,-We2^+ = Exp [-i2Tr/p Jei.ea] Wg^.+ Se^ - 

(Se,.)''= (We,.)'' = i (2) 



We require four operators on each edge, rather than two, 
in order to be able to associate our spin variables with 
fiuxes: in what follows, we may view S+ as measuring 
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S,,=2 



Se = ly 



FIG. 2. (Color online): This figure gives a selection of the con- 
figurations satisfying the vertex condition X]eG*v ('^eSe) = 
in a p = 3 Potts gauge theory. Configurations are represented 
by making edges with Se = 1 blue, Se = 2 red, and those with 
Se = Q dotted. The orientation of each edge is indicated by 
an arrow: on edges where the arrow is entering (leaving) the 
vertex, we take v = —1 (i' = 1). 



the electric flux in the x, y, or z directions, and S_ as 
measuring flux in the —x, — y, or —z directions; similarly 
We^± measures e**-'^^ with dl oriented in the 
or z direction. 

The Hamiltonian we will study has the form: 



(3) 



The three operators are defined as follows. The edge 
operator h is: 



^ m=0 



(4) 



he is diagonal in the spin basis, with eigenvalue I if Sg = 
0, and otherwise. This term therefore assigns an energy 
penalty of F to all edges with non-zero spin. 

The operator Bp acts on a plaquette P according to 

p-i 



P — J_ »• e,i/^ 
eedP 



(5) 



Here dP is the set of edges in the boundary of P, ori- 
ented clockwise relative to the x,y, or z normal, and i^^ 
is positive for edges in dP oriented in the x,y, or z di- 
rections, and negative for edges oriented in the —x, —y, 
or — z directions (see Fig. fTl). The relations (p|) imply 

that ($p) simultaneously raises (lowers) the spin on 



every edge of P traversed in a positive (negative) direc- 
tion by m units. It follows that {^P^ = li a-nd the 

flux (or eigenvalue of $p) through the plaquette has p 
distinct values, e^"'^-^/^', 0p = Q...p- 1. Moreover, Bp 
assigns an energy penalty of 1 to states with nonzero flux 
<j)p 7^ 0, thus favoring plaquettes with zero flux. 
Finally, the vertex operator is 



m=0 \ee*V 



(6) 



where *V is the set of all edges entering the vertex V . 
Here v is negative for edges whose separation from V is 
in the — i, — or — z directions, and positive otherwise. 
Q assigns an energy penalty of 1 for each vertex at which 
Eee*y(''eSe) ^ (mod p). 

The Hamiltonian ([3| has a few key features worth em- 
phasizing. First, 



Qv.Bp 



Qv,K 



Qv,Qi 







(7) 



and Qv is conserved separately at each vertex. We may 
therefore restrict our Hilbert space to states for which 
Qv — 1, which is equivalent to studying the pure gauge 
theory analyzed in Ref.'s [53] and [^Hl This imposes a 
constraint on the values of Se about a vertex, illustrated 
in Fig. ([2|. For p — 2, for example, edges with Se = 1 
must form closed loops. If we include configurations that 
violate this constraint, the model ^ is equivalent to a 
gauge theory with Potts-spin matter sources that cost 
energy M at each vertex. 

Second, the Hamiltonian ^ has two solvable points. 
If we take A = 0, then the Hamiltonian is diagonal in the 
spin basis, with he assigning an energy cost to any edge 
with non-trivial spin. In this limit the ground state is 
the product state with all Se = 0; loops with Se = m are 
low-lying excitations, whose energy scales linearly with 
the loop length. 

The model can also be solved exactly for F = 0, since 



Bp, Bp, = Qv,Qv' =0 



(8) 



and the Hamiltonian consists only of commuting opera- 
tors. Here the ground state satisfies the condition <I>p = 
on every plaquette; low-lying excitations are plaquettes 
for which (f)p — 1,2, ...p — 1. Because 13 p raises the edge 
spins, for F = each edge has an equal probability of be- 
ing in any one of the p different spin eigenstates. These 
ground states, which are superpositions over all spin con- 
figurations for which Qv has eigenvalue 1 at every vertex, 
are one example of what is generally called a loop gas for 
p = 2, or a string-net for p > 2. 

Between these two solvable points, the model under- 
goes a single phase transition, in which the non-zero spins 
are confined by the proliferation of vortex loop defects. 
To understand the two phases, let us begin with the F = 



4 



exactly solvable point. Here the ground state |4'o)r=o 
satisfies: 

-Bp|*o)r=o = |*o)r=o i-c $p|'J'o>r=o |*o)r=o • (9) 

Using the relations in Eq. (|2|, it is easy to show that hg 
and Bp fail to commute whenever e € dP. Indeed, 

(10) 

This implies that acting on |5'o)r=o with (Se)™ creates 
a small vortex loop of flux —2Trmi>e/p around the edge 
e, thus exciting all plaquettes bordered by e. In the de- 
confined phase these vortex loops, though present micro- 
scopically in the ground state for F ^ 0, remain short. 
As r/A grows, eventually these vortex loops proliferate, 
and the system undergoes a transition into the confined 
phase. 

Because the Zp°**'' models are self-dual and undergo a 
single transition as a function of T/X,^^ one can formulate 
a similar, dual description of the same transition starting 
from the A = exactly solvable point. Here the ground 
state has = 0, satisfying: 

he|*o>A=0 = |*o)a=0 i.e. Se,+ |*o>A=0 |*0>A=0 • 

Acting on the ground state for A = with (^p)™ creates 
an excited state in which a small loop of edges bordering 
the plaquette P carry non-zero Potts spin m. In the 
confined phase, these loops (or nets, for p > 2) of spins 
remain small. As F/A shrinks, eventually these loops (or 
nets) proliferate, and the system undergoes a transition 
into the deconfined phase. 

Thus the phase transition separating the confined and 
deconfined phases can be viewed either as a proliferation 
of vortex loops (for F/A increasing) or a proliferation of 
non-zero spins (for F/A decreasing). This transition is an 
example of a transition that is outside of the scope of the 
standard Landau- Ginzburg paradigm: there is no local 
order parameter that takes on a non-zero expectation 
value as we cross the transitiorli^. Rather, we may probe 
the onset of confinement by the (non-local) Wilson loop 
operatoil^. To understand this order parameter, we must 
consider states in which Qy^ I^E") = at a pair of vertices 
Vi,V2- In the gauge theory language, we can think of 
these states as having a pair of test charges at Vi and V2, 
which we can create by acting with 

w(^^= n wL^.., (12) 

efceCi2 

where C12 is a curve connecting vertices Vi and 
(Here = 1 if C12 crosses the fc*'* edge in the x,y, or 
z direction, and —1 otherwise; the p — 1 possible choices 
of q correspond to the possible charges 1, ...p — 1 of the 
discrete gauge theory). In the limit F = 0, the energy for 
creating a pair of charges is always 2M, irrespective of 
their separation; throughout the deconfined phase their 



energy scales like |ri2|° at large separations. Conversely, 
as F/A —> 00 the energy of separating the charges clearly 
grows linearly with |ri2|; this scaling holds at large sepa- 
rations throughout the condensed phase. In other words, 
the Potts gauge theory has p distinct charges (includ- 
ing 0), which are deconfined in the uncondensed phase, 
and confined in the condensed phase. The Wilson loop 
operator diagnoses this change in energy cost. 

For our purposes a closely related diagnostic - the 
topological orders of the two phases - will prove more 
practical. For A = 0,F > 0, the ground state is the 
unique spin-polarised state with Sg = 0, in which hg 
has eigenvalue 1 on each edge. In contrast, we will soon 
see that when A > 0,F = the ground state has a de- 
generacy of p^ for periodic boundary conditions. This 
change in the ground state degeneracy indicates that the 
uncondensed phase is topologically ordered, while the 
condensed phase (in which edge spins are confined) is 
not. Topological order cannot change without a phase 
transition's]^ so the system must undergo a phase transi- 
tion as F/A increases from to 00. 

The phase diagram of the Zp models has been studied 
in detail both numericalljESl and through large p series 
expansion^^Hm These results confirm that there are in- 
deed two phases, each of which is adiabatically connected 
to one of the solvable points discussed above; for p > 2, 
in 3D these are separated by a first-order transition. 

It is worth noting that in 2D these models are dual 
to the more familiar transverse-field Potts modeF^. In 
this dual description the phase with large F/A is param- 
agnetic, with spins tending to align with the transverse 
field. Large F/A corresponds to the ferromagnetic phase. 
Edges with Se ^ correspond to domain walls in the 
Potts ferromagnet. This relationship, and its relevance 
to the 2D cousins of the topological lattice models we will 
treat in the next section, are discussed in detail in Ref. 

m 



B. Topological order in the Potts gauge theory 

Hamiltonians exhibiting different topological or- 
ders necessarily represent distinct phases of matteif^Sl. 
Throughout this work we will shed light on the phase 
diagrams of various systems by distinguishing the topo- 
logical orders at exactly solvable points in the phase di- 
agrams. Our primary test of topological order is the 
ground state degeneracy on non-simply connected mani- 
folds, which we now calculate for the Potts gauge theory. 

Let begin in a ground state at the exactly solvable 
point F = 0, and consider the effect of the operator 

Wc, ^ n We,+ (13) 

where is a closed curve encircling the system in the 
fi = x,y,z direction (see Fig. For F = 0, W^^ com- 
mutes with the Hamiltonian. Thus for any ground state. 
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The topological ground state degeneracy is n^^^, where 
Ufree is the number of deconfined charges, and Nc the 
number of non-contractible curves. Indeed, the operator 



where is given in Eq. (12), and the curve from Vi 



Wc^ defined in Eq. (|13[) can be understood as wj^\ 



to Vi winds around the system m the z direction. 

What is the fate of this ground state degeneracy as the 
system crosses the phase transition? As A — >■ 0, config- 
urations which have any overlap with the original non- 
trivial ground states ^ 1 (for a = x,y ot z) become 
highly excited states with energy TLz, which diverges in 
the thermodynamic limit. This is because any such con- 
figuration has a non-contractible cycle (one which wind 
around the system in the x, y, or z direction) on which 
Se 7^ 1- Thus as soon as the phase boundary is crossed, 
these configurations are eliminated from the low-energy 
Hilbert space, and the model has no topological order. 



FIG. 3. (Color online): Ground state degeneracy in the pres- 
ence of periodic boundary conditions: With periodic bound- 
ary conditions, there exist operators Wc^, that raise spins 
on all edges along a non-contractible curve , and commute 
with the Hamiltonian. (The path Cz is shown in figure.) Act- 
ing with Wo,, changes the ground-state sector. The operators 
, which measure all the edges cutting a plane perpendic- 
ular to the fl direction, comprise a complete set of labels for 
the resulting ground states. In the figure, jj. = z and the edges 
between z = and z — 1 are measured. 



we must have: 



WcJ*o) 



(14) 



with also a ground state. We will find that the sys- 
tem is topologically ordered because these ground states 
are physically distinct. 

To see this, consider the operator 



(15) 



eec^ 



where Cz contains all edges connecting the planes z = 
and z = 1. commutes with $p for every P: if the 

latter raises the spin on one edge in C, it simultaneously 
lowers it on an adjacent edge, preserving the value of the 
product. However, 



(16) 



as there is exactly on edge on which both operators 
act. The states (WcJ |*o), m = 0,...p- 1 there- 
fore all have distinct eigenvalues under v^j (and similarly 
for X and y). We will call the ground state for which 
,V(j , all have eigenvalue 1 the trivial ground 
state. 



C. Relation to Zp lattice gauge theory 

It is worth noting that for p > 3 the operators ([4]) and 
([5]) differ slightly from their counterparts in the conven- 
tional Zp gauge theory, which are given by: 



h, = ^ S, 



^(^ 



■) 



(17) 



$1 



where $p 



Far from the transition, this model has the same phys- 
ical properties as the Zp°'*'' model: for F = 0, the ground 
state is flux-free (0p = 0), and at small F/A there is a 
phase with p deconfined "charge"-type excitations and 
a topological ground-state degeneracy of p'^ in periodic 
boundary conditions. For F/A very large, similarly, the 
ground state is a product state Se = 0, in which all spin 
labels are confined and there is no topological order. 

The main difference between these models and the 
gPotts niodels we focus on here is at intermediate F/A: 
for p > 5, in addition to the two phases of the Zp"***^, the 
Zp gauge theory has a third phase with gapless gauge 
excitations; this becomes the Coulomb phase of electro- 
magnetism in the limit p — o(^^^. These early analyses 
further suggest that the two transitions into this gapless 
phase are second order .fSSl 



III. CONDENSATION IN 3D CONFINED 
ABELIAN WALKER- WANG MODELS 

In this section, we will discuss a transition similar to 
the confining transition of the Zp models of the previous 
section in a very different family of 3D topological spin 
models: the Confined Abelian Walker- Wang (CAWW ) 
models discussed in Ref.'s [TS] and [2H Before entering 
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FIG. 4. This figure illustrates how the 6-valent vertices on 
a cubic lattice are split into a collection of trivalent vertices. 
The spin degrees of freedom living on each edge are repre- 
sented by black dots. 



into the details, it is useful to compare CAWW models to 
those found in the previous section. Like the Potts gauge 
theories of the previous section, the CAWW models have 
two phases separated by a first-order transition, which 
can be viewed as a condensation of vortex loops. In the 
uncondensed phase the models have loop gas (for p = 2) 
or string-net (for p > 2) type ground states, with decon- 
fined edge spins; in the condensed phase the non-trivial 
spins are confined, and (in the limit F — > oo) ground 
states are product states with = on every edge. 

Unlike in the previous section, we will find that nei- 
ther phase has conventional topological order. In order 
to distinguish them using criteria such as Wilson loops 
or ground-state degeneracy, we will have to study the 
models in the presence of boundaries. This is a conse- 
quence of the fundamentally different type of topological 
order found in the Walker- Wang models, which we will 
call surface topological order. 



A. Walker- Wang Hamiltonians 

We begin with a brief review of the Walker- Wang 
Hamiltonians discussed in depth in Refs [T5] and In 
this section we will focus on a subset of these, which we 
call Confined Abelian Walker- Wang (CAWW ) models as 
they are closely related to the Abelian lattice gauge the- 
ories discussed in the previous section, but do not have 
any deconfined bulk excitations (even if we allow vertex 
defects). For technical reasons the models are most eas- 
ily discussed on a lattice for which all of the vertices are 
trivalent; we thus deform the cubic lattice by separating 
each of its hexavalent vertices into three trivalent ones, 
as shown in Fig. [4j This deformation has little impact 
on the physics, so readers may safely imagine that the 



lattice is cubic unless otherwise notified. 

Like the Potts gauge models of Sect. |lT) the CAWW 
models begin with a Potts spin s = 0, 1, ...p on each edge 
of a 3D lattice. Schematically, the Hamiltonian has the 
same form as ([3|: 



H 



with 



and 



Bp, Bp, = Qv,Qv' = Bp,Qv =0 



hp, hp 



0. 



(18) 



(19) 



(20) 



This guarantees CAWW models share two key features 
with Potts gauge theories. First, as Qv commutes with 
all operators in the Hamiltonian, we may consider the 
limit M — oo, in which the constraint Qv = 1 is im- 
posed at each vertex. Second, the model can be solved 
exactly both for T/X — and for F/A = oo. Using these 
two limits, we will be able to characterize exactly the 
topological order of the two phases of this system. 

For the CAWW Hamiltonian, we take Qv to be given 
by Eq. ([Hi, and he to be given by Eq. Q. The only 
difference between these models and those of Section |ll] 
is in the definition of the plaquette operator. For the 
Walker- Wang models 



It. 



(21) 



where $p is the operator that raises all spins on the edges 
of P, exactly as in Eq. (jsj. The difference between Bp 
as given in ([s]) and (21 ) is the presence of the additional 
operator Qp^m- This operator is diagonal in the spin 



basis, with possible eigenvalues e 



iTzn/p 



0,...2p- 1; 



its value depends only the values of the spin variables on 
certain edges of P, and those on the 'legs' of the plaquette 
(i.e. edges that share a vertex with a pair of edges of P). 

The precise form of the operator Op.m for general p is 
not central to our discussion, but it is given in Appendix 
[K\ For our purposes, it suffices to note that, as shown in 
Appendix lAl 



e 



P,ni 



$p e 



p.i 



= 1 



= $ 



e 



P,m 



(22) 



The spectrum of Bp is therefore identical to that of 
the plaquette term in the 2^°'*^ gauge theory: Since 

0p,i) ~ 1' the eigenvalues of are 
A = e*'*^', (/)p = 27r n/p, < n < p. Further, 

1 



Bp = 



P 



E 



4>p e 



(23) 
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FIG. 5. (Color online): The plaquette operator for the p = 2 
CAWW Hamiltonian, illustrated for the three different pla- 
quette orientations of the point-split cubic lattice. The choice 
of red and blue edges used in Eq. ( 28 1 is indicated by the 



colors (and the letters R and B). All edges shown in the fig- 
ure enter into the plaquette operator: edges in the plaquette 
boundary dP are shown in bold; the remaining 10 "external" 
edges are in *P. 



Bp therefore has eigenvalue 1 on states where (f>p — 
(which we can interpret as states with trivial magnetic 
flux), and eigenvalue on states with (f>p ^ (which are 
thus states carrying non-trivial flux), exactly as in the 
Potts gauge theory. 



1. Example: p = 2 

For concreteness, consider the example oip — 2. This 
describes a system with an Ising spin — ±1 on each 
edge. We may represent 



iv 



n 



States that satisfy the condition 



(24) 



(25) 



are thus states with down spins (in the basis) on either 
or 2 of the three edges entering each vertex. In other 
words, for Af — ^ oo, we keep only configurations for which 
down spins form closed loops on the lattice. 

The plaquette operator for the Z2°"^ (or Ising gauge) 
Hamiltonian is 



B 



Potts 



1 



eedP 



(26) 



Since the non-identity term acts by simultaneously flip- 
ping all of the spins around a plaquette, the product 
over all plaquettes Y\p Bp contains matrix elements (of 
equal amplitude) connecting every possible configuration 
of closed loops with every other. (Here we are excluding 



loops encircling non-contractible curves, such as those 
mentioned in Eq. (13 1, from our definition of "closed" 



loops). The ground states satisfy 



Bp\^) 



(27) 



and are therefore equal-amplitude superpositions over all 
possible closed loop configurations on the lattice, with 
relative coefficient 1. Excited states can be formed by 
acting with Se,+ — cr| on some number of edges, changing 
the sign of the coefficients of configurations in which (t| = 
-1. 

The plaquette operator for the Walker- Wang Hamilto- 



Bp 




■Sj.-\-S 



(28) 

where *P is the set of all "legs" of the plaquette P, and 
r, r', are edges in the boundary of P, as shown in 
Fig. (|5|. Here Se — 5 (1 ~ "'I) = 0, 1 is the value of 
the PotIs spin on edge e. It is not hard to check that 
(B^W)2 ^ ^ww^ provided that Eq. ^ is satisfied at 

every vertex. This implies that B^^ is a projector with 
eigenvalues and 1. 

The non-identity component of acts by simul- 

taneously flipping all of the spins around a plaquette, 
and multiplying the resulting configuration by a phase 
of ±l,±i which depends on the spin configuration in 
*P and dP. The ground state is thus a superposition 
over all possible configurations of closed loops, with co- 
efficients equal in amplitude but with different relative 
phases ±l,±i. As for the Zf"**'' model, we can create 
plaquette excitations by changing the relative phases of 
these loop configurations. 



B. Topological order in the Confined Abelian 
Walker- Wang models 

We begin by reviewing the defining features of the 
CAWW models for F = 0. Since B^^ commutes with 
Qy, in this limit the models are exactly solvable. The 
resulting spectrum is discussed rigorously in Ref. [22| our 
discussion here is more qualitative, and we refer readers 
to that work for a detailed treatment. 

The presence of <d in the plaquette term of the CAWW 
model (Eq. (21 1 ) fundamentally changes the nature of 
the low-lying excitations, as well as the model's topologi- 
cal order. Unlike the 2^°'**^ gauge theories, these models 
have no deconfined excitations in the bulk, even in phases 
where the spin labels are deconfined and the ground state 
is a string-net. Commensurate with this, these 3D mod- 
els have a non-degenerate ground state, even for periodic 
boundary conditions, as well as a vanishing topological 
entanglement entropy. However, on a lattice with bound- 
aries, the models do have deconfined point particles on 
the surface and a ground state degeneracy sensitive to 
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the topology of the surface. Hence, the models do exhibit 
topological order, but it is associated with their surfaces, 
rather than with the bulk. 

It is worth noting that in CAWW models, generically 
the surfaces are described by topological orders that can 
be realized by a purely 2D system, suggesting that the 
surface topological order can be destroyed without un- 
dergoing a bulk phase transition. In the absence of extra 
symmetries, therefore, in most cases we know of no crite- 
rion which can differentiate the CAWW model from the 
topologically trivial confined phase. We will return to 



this point in Sect. VII 



Let us briefly review why the bulk ground states and 
spectrum differ from those of the Z^"*'** models. Essen- 
tially, the CAWW models have a unique ground state 
in the presence of periodic boundary conditions be cause 



there is no analogue of the operator Wc (Eq. 13 1 that 



commutes with B'^^^^ for every plaquette P and toggles 
between the different topological sectors. 

The origin of the differences in the low-lying excita- 
tions between the CAWW and Zp"**** models is similar. 
For r = 0, Eq. ( 18 ) indicates that excitations in the 



CAWW model can violate either the vertex term (as for 
the charge excitations in the model), the plaque- 

tte term (as for vortex loops in the Z^"**'' model), or 
both. Plaquette violations behave essentially identically 
in both models: they must form closed vortex loops, with 
an energy cost proportional to the loop length. In the 
CAWW models, however, it is not possible to create a 
pure vertex violation (see Appendix [A|). Rather, the low- 
lying bulk excitations in the CAWW models consist of 
vortex loops (aka closed loops of plaquette defects), and 
open flux tubes with electric charges at their end-points 
(literally, a pair of vertex violations separated by a line 
of plaquette defects). 

For instance, in the p = 2 CAWW discussed 
above, an operator that raises the spins on two legs 
of the plaquette P changes the eigenvalue of 9 



changed by - 



(Eq. 



28 1. (Specifically, it is 



-1 or ±i, depending on whether any of the 
r,r',b or b' edges are also flipped). Thus any operator 
that only raises the spins on these edges fails to commute 
with Bp. More generally, we might consider a modi- 
fied operator, which both raises the spins and multiplies 
the result by a configuration-dependent phase. It can be 
shown*^ that no such operator commutes with Bp for 
all P (with periodic boundary conditions). This means 
both that there is no operator that commutes with 
the Hamiltonian, and that vertex defects are confined. 
The latter follows because creating a pair of vertex de- 
fects at vertices vi and V2 requires raising the spin by 1 
along a continuous set of edges connecting vi to V2, and 
this raising operator necessarily leaves a string of violated 
plaquettes in its wake. In other words, the operator that 
creates a pair of "charges" (vertex violations) also creates 
a string of plaquette violations separating these charges, 
which we interpret as a magnetic flux tube. 



The previous discussion is relevant to systems with (for 
example) periodic boundary conditions. In systems with 
open boundaries, the CAWW models have additional ex- 
citations that are localized at the surfaces, and a ground- 
state degeneracy that depends on the surface topology. 
A useful intuitive picture of this is that the flux tube 
connecting the pair of charges can sit just above the sur- 
face, where it does not cost energy. (In technical terms, 
it is possible to define an operator which creates a pair of 
vertex defects on the surface without violating any pla- 
quettes. This operator acts by a combination of raising 
the spin along a series of edges on the surface, and multi- 
plying by a configuration-dependent phase.) These sur- 
face "charges" in fact behave like 2D charge-flux bound 
states, in the sense that they their mutual statistics are 
those of a J' = 1/p Laughlin state. (In more technical 
terms, these surface states have the topological order of 
an Abelian Chern-Simons theory). 

Since the charges are deconfined, a process that cre- 
ates a pair of charges, moves them along the surface, and 
re-annihilates them does not change the system's energy. 
It follows that if there are non-contractible curves on the 
surface^ there is an operator analogous to Wc that raises 
spins along a closed non-contractible loop and commutes 
with the Hamiltonian. In this case, however, such oper- 
ators do not necessarily commute with each other. For 
example, if we impose periodic boundary conditions in 
z only, the boundary of the system is a torus with two 
non-contractible curves, and we may define the opera- 
tors Wc^ , Wc^ . However, these two operators do not 
commut^^, and thus are not simultaneously diagonaliz- 
able. The ground state degeneracy is therefore p. Again, 
this matches the ground state degeneracy of the v = l/p 
Bosonic Laughlin state. 



Condensation transitions in the Confined 
Abelian Walker- Wang models 



We are now ready to consider the Hamiltonian ( 18 ) for 



r > 0. In particular, we wish to understand the structure 
of its phase diagram, and the nature of any phase tran- 
sitions. Much of this information will be inferred from 
our knowledge of the phase diagram of the Zp model 
([3]) -despite the significant differences between them for 
r = 0, the phase portraits of the two models are very 
similar. 

First, note that the effect of the transverse field term 
he is the same in both models: it penalizes any edge 
with Sg 7^ 0. In the limit F/A — ?■ 00, the ground states 
of both models are simply the product state Sg = 0; for 
large F the non-zero spins (and therefore all charges) are 
confined, and cannot form extended loops or nets. Thus 
for sufficiently large F/A the two models are in the same 
(topologically trivial) phase. 

It is worth pointing out that the mechanism for con- 
finement in the CAWW models with F = is fundamen- 
tally different from the mechanism for confinement in this 
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trivial phase. To see this, let us consider the p = 2 model 
of Sect. |III A 1 In the trivial phase, the transverse field 
term '^^<J^ penalizes any edge with = —I, meaning 
that only short loops of these spins can occur. In the 
Walker- Wang phase with F = 0, conversely, there is no 
penalty for an edge with — —1, and indeed the ground 
state of this model is a superposition over loops of all 
sizes. Hence if we view edges with cr^ = — 1 as occupied 
by non-trivial electric flux, electric flux is deconfined in 
the uncondensed phase, and confined in the condensed 
phase. However, there are no deconfined point-like exci- 
tations in the uncondensed phase, as the model does not 
admit any excitations that carry only electric charge. 

Second, the transition into this trivial phase is identi- 
cal in both models. As we show in detail in Appendix 
[Cj for the edge and plaquette degrees of freedom that 
are not conserved by (18), all correlation functions must 



be the same. Essentially, this stems from the fact that 
the commutation relations between the non-commuting 
operators in ([3| and ( 18 ) are identical. Specifically, for 
eedP, 



(*pep,i) 



(29) 



so that in both models y^ej raises the flux eigenvalue 

of (^^pQp^ij by 27rn/p, creating one of n distinct types 

of vortex loop excitations. This can be used to show that 
the amplitude for finding a given configuration of vortex 
loops in the ground state at a particular value of F/A 
is thus equal in both models (see Appendix [C]) , which 
ensures that the transitions are the same. 

An alternative way to understand this result is the fol- 
lowing. The only difference between the two models is 
the presence of Op,m in the plaquette projector, which 
changes the relative phases of different spin configura- 
tions in the ground state. However, these phases do not 
affect the expectation values of operators diagonal in the 
spin basis. (We take different spin configurations to be 
orthogonal). Rather, they are only sensitive to the prob- 
ability that Se takes on some particular value- which is 
the same in both models (provided we work in the triv- 
ial ground state sector). Put another way, the physical 
difference between the two models is in the charge sec- 
tor, with charges confined to the ends of flux tubes in 
the CAWW case. The transition, on the other hand, is 
in the purely magnetic sector, involving a condensation 
of vortex loops. From this perspective, it is unsurprising 
that the nature of the phase transition is the same in 
both models. 

As in the previous section, for p > 3 we can contem- 
plate a variant of our Hamiltonian, with hg replaced by 
( 17), and the plaquette operator given by ^pQp^i + h.c. 



Because the arguments given above apply equally well 
here (see Appendix [C|), the thermodynamic properties of 
the two models are the same. Hence for p > 5 in the 
"clock" variant of Eq. ( jTs] ), the two confined phases are 
separated by a critical. Coulomb-like phase. In this phase 



neither the vortex loops nor loops of non-trivial spin need 
have a finite energy cost - suggesting that the charge (or 
vertex-violating) excitations also become deconfined in 
this regime. 

Before discussing more general transitions in these sys- 
tems, it is worth pausing to consider what the results 
of this section tell us about the phase portrait of the 
Hamiltonian (18). In the Potts gauge model, the phase 



transition separates systems with different topological or- 
ders, which are therefore necessarily distinct phases. In 
the Walker- Wang system, however, the phase transition 
connects two systems with trivial bulk topological or- 
der. This raises the question of whether the confined and 
deconfined "phases" in the CAWW are truly distinct - 
or whether, like the liquid-gas transition, this first-order 
transition is an artifact of the particular trajectory that 
our models take through the phase diagram. 

For concreteness, let us focus on the case p = 2. In 
addition to topological order, we might identify two dis- 
tinct phases by identifying a symmetry that is broken 
by the CAWW phase but not by the trivial phase, or 
vice versa. However, this model possesses only lattice 
symmetries and the discrete symmetry PT, which sends 
X — ^ —X, y — > —y, z — > —z and complex conjugates the 
Hamiltonian. Because of the nature of the plaquette op- 
erator, both PT and the lattice symmetries are symme- 
tries of the system for all values of F, A, while P and T 
are individually broken unless A = 0. 

Thus for p = 2 the CAWW model cannot be distin- 
guished from the trivial confined phase either on the ba- 
sis of symmetry or on the basis of topological order. This 
strongly suggests that there is a deformation of (18) al- 
lowing us to connect the two limits F = and A = 
without undergoing a phase transition. The arguments 
in Appendix [C] ensure that such a deformation must in- 
clude terms that create vertex violations, as otherwise 
the free energy is identical to that of a model whose pla- 
quette term contains only the operators $p, in which 
the two limits have different topological orders. Further, 
as the phase transition is first order, it will persist for 
arbitrary small deformations away from ( 18 ) - possibly 



ending at a second-order point at some finite value of the 
perturbation strength. 



D. Diagnosing vortex condensation 

How do we discern the behaviour of the CAWW (F = 
0) from that of the condensed phase (F = oo)? The 
condensation transition does not appear to involve the 
symmetry breaking of an order parameter and, while the 
presence of bulk topological order distinguished the Potts 
gauge theories from the trivial paramagnet, we cannot 
say the same of CAWW models which have no bulk topo- 
logical order. One obvious difference between the phases 
is in the presence of many vortices, which can be detected 
by the operator 
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^5 = n ^p^p-.i (30) 

Pes 

Os is formed by acting with the plaquette operators on 
open surface 5, with area A and perimeter L. Near 
the r = exactly solvable point, the operator obeys 
a perimeter law {Os) ^ e~^^ because vortex loops are 
short, so can only intersect the surface if they lie near its 
boundary. On the other hand, near the F — !■ oo point, 
{Os) ~ e"'^'^, because vortex loops are larger and pro- 
liferate freely. While the change from area to perimeter 
law does not imply that the F = and F = oo points 
constitute distinct phases of matter, it does allow us to 
quantitatively distinguish the two regimes. 

In the Potts gauge theory the fundamental Wilson loop 
along dS, which can be written as O^'^^^^ = Hpes ^p^ 
plays the same role as Os in tracking vortex condensa- 
tion. However, if used in the CAWW case, this operator 
will exhibit area law behaviour in both the un-condensed 
and condensed regimes. This is because, even at F = 0, 
the magnetic flux measurement $p fluctuates wildly - a 
defining feature of CAWW models is that only the spe- 
cial combination ^pQpi of electric and magnetic flux 
measurements is stationary on the ground state. 



IV. OTHER TRANSITIONS IN ABELIAN 
MODELS 

Thus far, we have identified a set of first order tran- 
sitions in the Z^"**'' and CAWW models which connect 
a given 3D string-net state to a trivial phase, in which 
all spin labels are confined. We now turn to a rather 
more interesting situation, in which only some of the 
spin labels are confined in the condensed phase. Here 
we will find a surprising difference between the Z^"**'' 
and CAWW models. In the Zp°**'' case, a transition that 
confines spin labels is a confining transition: the number 
of deconfined charges in the gauge theory decreases, as 
does the ground-state degeneracy; at long wavelengths 
the confined phase may be described by an Abelian lat- 
tice gauge theory with a smaller gauge group than the 
deconfined phase. For the CAWW models, however, in 
the uncondensed phase there are no deconfined charges, 
and we might guess that the transitions are necessarily 
between two phases with trivial bulk topological order 
and no deconfined point-like excitations. We will find - 
somewhat counter-intuitively - that this guess is incor- 
rect: a phase transition which confines a subset of the 
spin labels in an CAWW model can lead to deconfine- 
ment in the charge sector. Correspondingly, we will find 
that the ground-state degeneracy grows as we cross the 
transition into the confined phase. We will give a physical 
interpretation of this in Sect. |IVB| 

How do we obtain transitions that confine only some 
of the labels? The key is that Eq. Q is not the only 



possible choice of he. More generally, if p is divisible by 
m we can take 

p/m—l 

h(™) ^^Y. (Se,+) (31) 

which has the effect of proliferating vortices with fluxes 
2tt/p X {m, 2m, . . . ,p—m}. The operator hi™-* has eigen- 
value 1 on edges where Se is a multiple of p/m, and 
otherwise, and therefore assigns an energy penalty F for 
every spin that is not a multiple of p/m. Thus even 
as F — > oo, some spins remain deconflned. The result 
for the ZP°"" models is a condensed phase that has the 
bulk topological order of a Zp/m gauge theory. For the 
CAWW models, depending on the values of p and m, we 
find that the condensed phase can have bulk topological 
order, surface topological order, or both. A table classi- 
fying the topological orders of these phases can be found 
at the end of Sect. HVBl 

As in the previous sections, "clock" -like versions of the 
Hamiltonians that we discuss also exist, and for m > 5 
these will exhibit gapless phases at intermediate values 
of F/A. Though this possibility is certainly intriguing, 
our focus here will be on the possible phases at large F, 
where the two models are in the same phase. 

For readers familiar with TQFTs, we note that not all 
Abelian Walker- Wang models are completely confining 
at F = 0. In Appendix [D| we classify the possible con- 
densed phases in arbitrary Abelian Walker- Wang models. 
Like the CAWW models discussed here, the structure of 
their phase diagrams is identical to that of an appropriate 
^Potts gauge theory. 



A. Variants on the 3D Potts gauge transition 

For pedagogical reasons, we will first describe these 
transitions in the Zp°'*^ models. If we replace he by 

hi™^ in Eq. |3| and study the resulting phase diagram, 
we again find two phases separated by a first-order phase 
transition. We will show momentarily that an effective 
Hamiltonian describing the transition can be mapped 
onto that of a system with p/m spin states and a trans- 
verse field hi^\ such that these transitions are always 
first order. However, the topological order of the phase 
at large F/A is different for every m: it is that of an m- 
state Potts gauge theory. Thus with a Hamiltonian of 
the form ^ and the modified definition of hi™'' , we can 
describe a family of transitions between our p-state Potts 
gauge theory and an m-state Potts gauge theory, for any 
m that divides p. 

We now explain why the phase transition is still of the 
type described in the previous section. Since 

Se,+) $P = e^2™/P$p (Se,+) (32) 
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we have 



p/m 



(33) 



Because 



p/rn 



commutes with ah other terms in the 



Hamiltonian, we may restrict our attention to states 
for which 



^ m— 1 
m 5Z \ 



np/m 



(34) 



To ensure that this condition is satisfied by aU ground 
states even for A = 0, it is convenient to modify the 
plaquette projector somewhat: 



XBf 



P 



p/m— 1 

A + A 

i=0 



^ 7n — l 



np/r, 



(35) 

The term proportional to A commutes with ah operators 
in the Hamihonian; thus we m ay consider the hmit A — >■ 
oo, where the condition (34 1 is always satisfied.-^ On 
states satisfying Eq. (34), we have: 



p/ni~l 



Bp\^) = - E ('^^) I*) 



(36) 



i=0 



Within this sub-space, then, the Hamiltonian for the 
transition is equivalent to a Hamiltonian of the form ([3| 
for a p/m state spin. 

Thus the energetics of the transition are clearly those 
of a transition with p/m-state spins and hi^''. The only 
remaining subtlety is whether the number of states at 
each energy is the same in both systems. We show in 
Appendix [B] that though the total number of spin config- 
urations at any energy is of course larger in the p-state 
system, this multiplicity merely alters the normalisation 
of the ground-state wave function, and plays no role in 
the thermodynamics of the phase transition (at any tem- 
perature). 

Next, let us verify that the topological order of the 
phase for large F/A is indeed that of an m-state Potts 
gauge theory. From Eq. ( 33 1 , we see that spins with Sg = 
n p/m carry no energy cost. It follows, from arguments 



similar to those given in Eq's ( 13 1 to ( 16 1 , that the system 
has an m^ ground-state degeneracy in periodic boundary 
conditions'7'. Similarly there are m deconfined charges in 
the system, which lie at the end-points of string operators 
corresponding to the m deconfined spins. 



B. Variants on the transition in Abelian 
Walker- Wang models 

We now consider transitions in the CAWW models 
where the edge term is given by h^"") (Eq. (ml) ), for 



m a divisor of p. Again, this choice assigns no energy 
cost to the p/m spins s = 0,m, 2to, ...p — m, while all 
other spins are eliminated from the low-energy Hilbert 
space as F/A — >■ oo. Since h(™) is insensitive to the rel- 
ative phases in the ground state wave functions, the ar- 



guments of Sect. IV A and Appendix [C] can be combined 
to show that the phase transition is equivalent to that of 
a p/m state Potts gauge theory restricted to its trivial 
ground-state sector - and therefore must be first order. 

Our main interest in this section, then, will be to de- 
scribe the topological order of the condensed phase. Since 
the topological order is a property of the phase, it is con- 
venient to study it by identifying a point in this phase at 
which the Hamiltonian is exactly solvable. To do so, as 
in the Z^"*'*' case (Eq. 35 ) . for m > 1 we will modify the 
Walker- Wang plaquette projector according to: 



XBp = 



-t m— I , 
1 ^-^ / - \ np/rn . 

A 2^ [<Ppj Qp.np/m 

p/m— 1 m — 1 . , , 

2^ 2^ ] ^P,i+np/m 

i=l n=0 



(37) 



. The term proportional to A commutes with h(™) ir- 
respective of the presence of O: 8 and h(™) necessarily 
commute, as both arc diagonal in the spin basis, and also 

/ - \np/m 

[h(™), I $p J ] = 0. Thus the model can be solved ex- 
actly in the hmit A = 0, F^-oo,A>0. 

For A = 0, the plaquette projector changes the values 
of the spins only by multiples of p/m. Taking F — oo 
ensures that the low-energy Hilbert space contains only 
the m spins s = 0,p/m, ...p — p/m. Thus by re-labelling 
the spins — > Q,p/m — > 1, ...p/m — 1 — )• m — 1, we may 
express the Walker- Wang Hamiltonian in this limit as 



ri=0 



n t 



(38) 

with Qy, $p the vertex and flux operators of the m- 
state Potts gauge theory, and 9 p,„ defined by the action 
of the operator Qp,np/m on the reduced Hilbert space 
containing only spins that are multiples of p/m. 

To proceed further, we must understand the basic 

properties of Op.n. In particular, we wish to determine 
whether there is an operator analogous to (Eq. ( 13 ) 
) that commutes with the Hamiltonian (38). If so, the 
system will have multiple ground states, distinguished by 
their different eigenvalues under the action of the opera- 
tor i>c in Eq. (15) (which does commute with _ffww)- If 
not, as shown in Ref. 1221 the ground state is unique. 

The eigenvalues of Op,„ can be deduced from those of 
Qp,np/rm which are phases depending on the spin con- 
figuration of the external edges of P, as well as the 
spin configurations of certain edges on P. In the low- 
energy Hilbert space containing only spins that are mul- 
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tiples of p/m, the possible eigenvalues of 6p.„p/,„ are 



„iiv(r /m){p/m) 



, where r = 0, 1, 



1. 



There are four cases to consider here. First, \ip/m £ 
2Z, then e*'^(''/™)(p/") ee 1 for ah r. In this case, one can 

show that 8 acts as the identity operator on all states 
obeying Qyl^*) = l^*). If p/m^ is an even integer, the 



Hamiltonian (38) therefore describes an m-state Potts 



gauge theory, with a topological ground state degeneracy 
of in periodic boundary conditions. 

Second, if p/m2 e 2Z + 1, then e*'^('-/™)(p/'") = (-l)'". 
In this case, we can still define an operator for each 
non-contractible C that commutes with the Hamiltonan 
for any value of A, such that the ground-state degen- 
eracy remains in periodic boundary conditions (see 



Appendix 
Appendix 



Al for a proof of this). In fact, we show in 
Elthat for M < CXI this model describes an 
TO-state Potts gauge theory coupled to fermionic matter 
sources. 

Third, if m is not a factor of p/m, then the only value 
of r G {0, ...TO — 1} for which Op^Wm has eigenvalue 1 is 
r — 0. As we show in Appendix [A] this means that the 
Hamiltonian (38) describes an TO-state CAWW model, 
with trivial topo ogical order. 

Finally, suppose pjrr? — ajb, with a and h relatively 

prime, and & > 1 a factor of to. In this case Q has b 
eigenvalues ±1 for r ^ tb, t — 0,1,. ..to/6 — 1. If 

has eigenvalue 1 for r 7^ 0, the eigenvalue of 6 is unaf- 
fected by replacing s — > s + r on any edge (see Appendix 

A). Therefore the operator ^M^c^ commutes with the 
"Hamiltonian. Less trivially, one can show that provided 
Qyl^*) — 1 4') (which is true in the ground states, pro- 
vided M is positive), ^W^c) also commutes with H for 

r corresponding to negative eigenvalues of Q. Thus in 
this case we find a model whose Hilbert space is de- 
scribed by TO-state Potts spins on each edge, with (in 
periodic boundary conditions) a ground-state degener- 
acy of {m/b)^ , where b divides to. The topological order 
of this phase is therefore different from that of both the 
TO-state Potts gauge and CAWW models. 

A summary of these possibilities, together with the 
topological order of the condensed phases, is given in 
Table IB 





Topological order of the condensed phase 


2n, n e Z 
2n + l, n£Z 
n/m (irreducible) 
n/b, b\m 


Zm gauge theory 
fermionic Z™, gauge theory 
Zm CAWW (modular) 
^mn/(26) Q^"^^ (non-modular) 



TABLE I. Topological orders of p state CAWW models in 
the condensed phase with transverse-field h''"-' . Our notation 
conventions are drawn from Ref. i38l 

What becomes of the excitations as we cross the 
boundary from an uncondensed phase to a condensed 



phase? Since all spins are now multiples of p/m, there 
are only to — 1 physically distinct vortex-loop creation 
operators 



(39) 



Writing the plaquette projector as in Eq. (37) allows 
us to separate the fluxes bound to our charges into two 
kinds: those whose energy cost remains finite (propor- 
tional to A) as A — >■ and those whose energy cost van- 
ishes for A = 0. The number of fluxes with zero energy 
cost in the confined phase is given by the number of spins 
r for which e^'^^ = ±1. Thus some particles be- 

come deconfined in the condensed phase, because they 
lie at the end-points of flux tubes whose energy cost van- 
ishes as A/F — 0. Physically, this is because these flux 
tubes become unobservable (and thus have vanishing en- 
ergy) once a subset of the spins have been eliminated 
from the theory. The deconfined particles are always 
fermions or bosons: the criteria e'^'' — ±1 ensures 
that the relevant charge-creation operators either com- 
mute or anti-commute. (For a more detailed discussion 
of the fermionic case, see Appendix [E|). 

It is also instructive to consider the fate of the decon- 
fined surface excitations in the condensed phase. There 
are three possibilities: a given type of surface anyon may 
become confined; may remain deconfined at the surface 
only; or may become deconfined in the bulk (as well as 
at the surface). 

To understand the first possibility, we need only know 
that the surface anyons exist at the end-points of a string 
of edges of appropriate spins. For example, in the p — 2 
example discussed above, a pair of surface anyons occurs 
at the end-points of a string of edges on which s = 1. In 
that case, in the condensed phase edges with s = 1 bear 
an energetic cost; hence the surface anyons are confined 
(in the same sense that charges become confined in the 
gPotts gauge theory) , and the transition is one to a state 
with trivial bulk and trivial surface topological order. 

If the anyon is to sit at the end-point of a string of 
edges with spin labels that remain deconfined, then that 
excitations of this type are deconfined on the surface in 
the condensed phase. Whether they are also deconfined 
in the bulk then depends on whether pjm} is an integer. 
If it is, then our "anyons" are in fact bosons or fermions, 
and in the condensed phase have trivial braiding with 
all of the deconfined surface anyons. If p/to^ is not an 
integer, then the condensed phase has excitations that 
are deconfined only at the boundary of the 3D lattice; at 
least one of these must be anyonic. 



C. Examples 

For concreteness, we will now consider one example 
from each of the four classes in Table (|l]). To obtain a 
physical picture of which excitations are deconfined in 
the bulk, it is helpful to classify the excitations in the 
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Walker- Wang models in terms of their "charge" (given by 
the eigenvalue of Qv at the violated vertices) and "flux" 
(given by the phase of the eigenvalue of $p.i6p.i on the 
line of excited plaquettes connecting the two vertices). 
In the uncondensed phase, a charge q lies at the end of 
a flux tube with flux (j) = 2'Kpq. The operator that raises 
edge spins of P by s essentially measures the Berry phase 
of an object of charge q = s/p around the flux through 
the plaquette P, penalizing states where this Berry phase 
is not a multiple of 27r. To understand why some charges 
become deconfined in the condensed phase, we observe 
that for A = 0, the plaquette term (37 1 contains only a 
subset of such possible measurements - those for charges 
q = np/m, n = 1, ...m — 1. In this limit, some flux tubes 
become physically unobservable: there are no longer any 
charges in the theory with which they have non-trivial 
Berry phase. This gives us a simple mnemonic for un- 
derstanding the topological order of the phase with large 

r/A. 

Armed with this intuition, we will consider each of the 
four classes (see Table |l]) in turn. For simplicity, we will 
discuss these in the solvable limit A = 0, F > 0, for which 
there is zero amplitude to create any edges with con- 
flned spin labels; however, the physical properties are 
not restricted to this special point but characterize the 
condensed phase. 

The first possibility is that the condensed phase can 
be described by a Z,„ Potts gauge theory. To see this, 
consider the example p = 8, m = 2. We can view the spin 
labels s = 0, ...7 as representing 7 possible electric fluxes 
that can exist on each edge of the lattice, associated with 
fractional charges qs = 0, 1/8, ...7/8. The possible mag- 
netic fluxes in the theory are then 27rn,n = 1,..7. The 
fundamental excitations in the uncondensed phase are 
closed flux tubes of flux 27rrt, and open flux tubes with 
flux 2TT{8q) which terminate at vertices with charge q. 

For A/F = 0, only the two spins s = 0,4 remain in 
the low-energy Hilbert space. If we keep only these two 
spins, the vertex condition requires that the number of 
edges with s = 4 at each vertex must be even - in other 
words, it reduces exactly to the vertex condition of the 
Z2 model. Further, if we keep only edge spins s = 0, 4, 
then 6)4 = 1, and the plaquette projector simply flips all 
spins around a plaquette. Thus our solvable Hamiltonian 
for the condensed phase is exactly that of Z2°*^^^ gauge 
theory. 

What has become of our charge-flux-tube bound 
states? In the condensed phase, the only deconfined 
charge is q — 1/2, which lies at the end of a flux tube 
of flux Stt. In the uncondensed phase, this magnetic flux 
tube is physically observable (and indeed costs energy 
per unit length), since it is measured by lines of elec- 
tric flux corresponding to charge q — 1/8,3/8,.... Elec- 
tric flux lines corresponding to q — 1/2, however, can- 
not distinguish between a flux 87r and a flux 0, and the 
(7 = 1/2 charges are the only deconflned in the condensed 
phase. The only other low-energy exctation in this phase 
is the closed vortex loop of flux 2{2n + 1)tt, about which 



a charge 1/2 has a Berry phase of tt. (All such vortex 
loops are physically indistinguishable in the condensed 
phase). This is exactly the spectrum of Z2 gauge theory. 

Next, consider p = 4, m = 2. The associated electric 
charges are q = 1/4, 1/2,3/4, bound to the ends of flux 
tubes of flux 27r(4g). For A/F = 0, the only remaining 
electric flux is 1/2, and the possible charged excitations 
are q = 1/2, bound to flux tubes 47r. These fluxes are 
physically undetectable by objects of charge 1/2, so that 
the charges are deconfined. The condensed phase also 
has vortex loop excitations, of flux 27r = Gtt. Indeed 
the only difference from the previous example is that the 
charge 1/2 has Berry phase 2tt, rather than Air, with its 
associated flux tube. As we show in Appendix |Ej this 
describes a Z2 gauge theory with fermionic sources. 

Third, consider p = 6,to = 2. Here again the electric 
fiux that remains for A/F = is 1/2; charge 1/2 excita- 
tions are now bound to the end-points of tubes of fiux Gtt. 
These fiux tubes have a Berry phase of tt with particles of 
charge 1/2, meaning that they are physically observable, 
and thus cost a finite energy per unit length. The charge- 
1/2 excitations (the only possible charges in this limit) 
are therefore confined. Indeed, this model is exactly the 
p = 2 CAWW model discussed in Sect. |IIIA1[ 

Finally, consider p — 8,m = 4. In this case there are 
four electric fiuxes that remain for A/F = 0, of strength 
1/4, 1/2, and 3/4. The corresponding charges are bound 
to fiux tubes of strength An, 8tt, and 127r, respectively. 
The first and last of these fluxes have Berry phase tt with 
objects of charge 1/4, and thus are physically observable. 
A fiux of 87r is physically unobservable, so that there is 
a deconfined excitation of charge 1/2 (corresponding to 
s = 4 in the original spin basis). According to our cri- 
teria above, this excitation is a boson. The excitations 
of charge 1/4,3/4 are confined. This model is therefore 
intermediate between the Potts gauge theories and the 
CAWW models: it has deconfined charge and non-trivial 
topological order, but the number of deconfined charges 
(and the ground state degeneracy) is less than the num- 
ber of deconfined spin labels. 



V. ISING GAUGE TRANSITIONS IN 
WALKER- WANG MODELS FOR SU(2)t, 

The CAWW models are only a subset of the possibili- 
ties that can be realized following the construction of Ref. 
1151 In this section and the next, we will discuss transi- 
tions in a different family of Walker- Wang models, whose 
low-lying excitations include non-abelian anyons confined 
to their surfaces. We will call these the SU(2)fc Walker- 
Wang models, as on lattices with boundaries, their sur- 
face states have the topological order of a chiral SU(2)fe 
Chern-Simons theory. 

In the SU(2)fc Walker- Wang models, each edge is en- 
dowed with a spin variable s = 0, 1/2, ...fc/2, where k 
is a parameter of the model. The vertex operator acts 
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according to 



where 



ifc \k 

l^ixj 



min{i-\-j.k—i—j} 

ix j = I 



(40) 



(41) 



Notice that if fc = oo, we recover the usual rules for the 
addition of spin angular momenta, and Qy has eigen- 
value 1 on states where the total angular momentum is 
conserved at the vertex V. Eq. (41) enforces a "de- 
formed" version of this condition, appropriate to models 
with a finite number of possible values of the total (de- 
formed) spin angular momentum on each edge.'^^ 

As in the Abelian models, the plaquette operator Bp 
is a superposition of operators that raise all spins on the 
edges around the plaquette by s (using the rules (41) 
), and simultaneously multiply the wave- function by a 
configuration-dependent complex coefRcient: 



k 

E 

2s=0 



(42) 



In this case we cannot separate the action into indepen- 
dent raising and phase operators, as in general raising all 
spins on the edges of P by s will create a superposition of 
configurations with different spin labels; each element of 
this superposition may have a different phase. Thus we 
denote by $p.s the combination of raising operator and 
phases. (The precise form of these operators is given in 
Refs. 



VIA) 



T5landl22l see also Sect. 
Following the procedure that we used in the Abelian 
models above, we will seek to identify possible edge terms 
(transverse fields) that we can add to the solvable Walker- 
Wang Hamiltonian to produce a phase transition. In 
the Abelian case, the transverse field operators he were 
sums of terms S" _|_ -I- S" _ , each of which has eigenvalues 
of the form cos (27rnse/p). For the SU(2)fe models, it 
is convenient to choose a slightly different form of the 
transverse field: 



hi™) I 



7r(2se + l)(2m+ 1) 
fc + 2 



(43) 



where Se,™ G {0, 1/2, 1, ...fc/2}, and n^^ ~ 
sin [7r(2se + l)/(fc + 2)]. This has a similar effect to 
adding a term of the form cos(47rjs/fc), but assigns 
slightly different energy penalties to the various edge 
spins. The advantage of the form (43) is that it obeys 
the rules (41) for combining spins: 



hi™) hi") I 



E 



(44) 



The excitations created by the operator hi™) is therefore 
a vortex loop "of spin m" , and can be combined according 
to the same rules (Eq. 41 1 we use to raise and lower the 
edge spins. 

For most choices of k,m, the transverse field will assign 
an energy cost to every non-0 spin, and the condensed 
phase is the "trivial phase" , in which all spin labels are 
confined. (If we choose m = 1/2, for example, this is true 
for all k). We will not discuss these transitions in detail, 
but note that they cannot be mapped onto the transitions 
discussed above: the commutator of hi™) with Bp is not 



a sum of terms of the form (29), and the condensing 



vortex loops do not behave like vortex loops in a Zp"**** 
(or related Abelian) model. 

There are, however, some choices of the transverse field 
term which do not lead to a transition into the trivial 
phase. These will be our primary interest here. In Sect. 
|VI| we will describe a general method to determine the 
topological order of the condensed phases for general to. 

We begin, however, with two examples. First, if we 
take TO = 1, 



hW|fc/2)-«^i|fc/2) hW|0) = ^i|0) 



(45) 



with Kl ~ sin 



37r 
k+2 



It is not hard to show that the 



eigenvalue of hi^) for all other spin states is smaller, so 
that this transverse field assigns an energy penalty to any 
edge label except 0, fc/2. In this case, we will have little 
to say about the behaviour near the transition. How- 
ever, we will show that the condensed phase is (1) the 
Z2 gauge theory if k is divisible by 4; (2) the Z2 gauge 
theory with fermionic charges if k = 4n + 2; (3) the p = 2 
CAWW model (Sect. |III A 1[ ) if fc is odd. The differ- 
ence between these cases arises from the different possible 
phases incurred by acting with $p_fc/2 on configurations 



with Se = 0, fc/2 only, as we discuss in Sect. VI 



A second interesting possiblity is to take to — fc/2. 
this case we have 



In 



|Se) = - COS {{2s, + l)Tr) = ("1)''= \Se) (46) 



such that — hi'^^^) assigns an energy penalty to any state 
where the edge e carries a half-integer spin. This case is 
the 3D analogue of the Ising transitions studied by Ref. 
|4]and[8l We will show that as in the 2D case, for every 
fc the phase transition is identical that of the Z2 gauge 
theory. The nature of the condensed phase, however, 
depends strongly on fc, as we shall see. 



A. Transitions in SU(2)2 

We begin by studying the simplest model in this family, 
the SU(2)2 Walker-Wang model. Here there are only 
three allowed spin labels s = 0, 1/2, 1, and since fc/2 = 1, 
the two choices (45) and (46) of transverse field coincide. 
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The Hamiltonian we will study is 



p 

7i 



(47) 



with the action of he given by Eq. ( 46 ) 
The action of the vertex term is: 



XXX 

gy|V^) = <(l <"■"■"> ("'1-11 U'i'l) (48) 



otherwise 



where it is understood that rotations of these combina- 
tions have the same eigenvalue under Qy. Since the 
eigenvalue of Qv is conserved at every vertex, we will 
restrict our attention to states where these eigenvalues 
are 1 everywhere. In this restricted Hilbert space, edges 
of spin-1/2 always occur in closed loops, while edges of 
spin-1 can either form closed loops, or open lines ending 
at vertices with two spin-1/2 edges. 

We have separated the action of the plaquette operator 
into two terms: 



(49) 



The term proportional to A changes the spin on each edge 
of the plaquette P by an integer amount, and thus com- 
mutes with the transverse field term (and in fact with all 
terms in the Hamiltonian). The term proportional to A 
changes the spin by a half-integer amount - and therefore 
changes the eigenvalue of he on all edges of P. (Walker 
and Wang'fli^ definition of the plaquette operator Bp is 
obtained by setting A = A.) 

In more detail, the action of the operators $p_s is as 
follows. $p^i raises the spin on every edge of P by 1, 
using the rules (see (41 1 ) 



1x0 = 1 



1x1 = 



1x1/2 = 1/2 (50) 



In other words, spin— 1/2 edges remain spin-1/2 edges, 
while edges of spin-0 and spin-1 are interchanged. (The 
operator also in general multiplies the wave- function 
by a non-zero configuration-dependent coefficient; we 
will have more to say about these coefficients presently). 
$P i/2 raises the spin on every edge of P by 1/2, using 
the rules (see (41 ) ) 



1/2x0 = 1/2 1/2x1 = 1/2 1/2x1/2 = 0+1(51) 

In other words, it turns edges of spin or 1 into edges 
of spin 1/2, and edges of spin 1/2 into a superposition of 
edges of spin and spin 1. Again, this action also multi- 
plies the result by a configuration dependent coefficient. 
Importantly, this coefficient is for any configurations 
where the eigenvalue of Qv has changed. This ensures 
that the plaquette and vertex terms commute. 



We now wish to prove the two assertions following Eqs. 
(45 1 and (46 1: first, that this transition is identical to the 



confining transition of the Z2 gauge theory, and second, 
that the condensed phase at large F/A can be described 
by the deconfined phase of Z2 gauge theory with fermionic 
sources. 

We begin with the second of these assertions. As in 
previous sections, to understand the behavior of the con- 
densed phase, it suffices to consider the Hamiltonian at 
A = 0,r > 0. In this case the spin on each edge is 
conserved modulo an integer, and each edge of spin-1/2 
carries an energy cost of F: spin-1/2 labels are confined, 
in the sense described in Sect. |llj Restricting our atten- 
tion to states containing only the deconfined spin-1 label, 
the effective Hamiltonian has the form 



H 



(52) 



where Qv penalizes states with an odd number of spin- 
1 edges entering the vertex V. In the restricted Hilbert 
space, $p_i simply interchanges all spin edges of P 
with spin 1 edges, and all spin 1 edges with spin edges, 
multipliying the result by a coefficient ±1. Apart from 
a possible sign in the action of the plaquette term, we 
recognize this as the Hamiltonian for a Z2 gauge theory, 
whose possible excitations are vortex loops of plaquettes 
on which $p^i = 0, and vertex violations (charges) where 
an odd number of 1 spins meet at a vertex. As we men- 
tioned in Sect. |IV B[ the fact that the matrix elements 
of the plaquette term are always real ensures that the 
charges are deconfined; the effect of the — signs in this 
case is to render the charges fermionic, as we prove in 
Appendix lEl 



The Hamiltonian ( 52 ) is, of course, valid only precisely 
at A = 0; however, the qualitative features of the Z2 
gauge theory describe the low-energy physics of the entire 
phase. Since spin-1/2 edges are confined throughout the 
large F/A phase, these loops can never appear at long 
length scales, and Eq. (52) is the effective Hamiltonian 



at long-wavelengths throughout this phase. 

Having understood the condensed phase, let us under- 
stand the nature of the transition. This is most easily 
done if we restrict our attention to states for which 



p,i 



IV) = IV) 



(53) 



i.e. states with no vertex violations, and no spin-1/2 
vortices. ($p,i does not measure the spin-1 vortices gen- 
erated by hi^''). In this subspace, for F = the model 
contains only one type of excitation: a plaquette defect 
visible only to $p.i/2, which costs energy A. This is pre- 
cisely the plaquette defect created by the transverse field 
term hi^\ which commutes with <l>p^i and anti-commutes 

with <&p4/2- Further, (he^') ^ 1> exactly as for the 
transverse field term in the Z2 gauge theory. Thus we 
might expect that the states in the Hilbert space defined 
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by Eq. (53) can be mapped onto states of the pure Z2 



gauge theory, by identifying the vortex loops created by 
hi^^ with the vortex loops in the gauge theory. 

It is instructive to construct such a mapping directly 
in the spin eigenbasis. On each edge, we may map the 
degrees of freedom of the SU(2)2 model to that of a Z2 
gauge theory via: 



— > 2se ( mod 2) 



(54) 



The reason that this is a sensible mapping is that at 
each vertex V, the number of half-integer spins entering 
V must be even; similarly in the Z2 gauge theory the 
number of edges with s — 1 must be even. Thus in 
the absence of vertex violations, the mapping is one from 
closed loops of spin 1/2 to closed loops of spin 1. Further, 

the transverse field term hi^^ = (— 1)^"" maps exactly to 
the transverse field term hg = (—1)^" of the Z2 gauge 
theory. Finally, the operator $p_i/2 interchanges integer 
and half-integer spins, while the plaquette term in the Z2 
gauge theory interchanges 8^ = and 1. 

The mapping (54), however, clearly fails to capture 
many of the features of the ground state of the SU(2)2 
Walker- Wang model. First, as a mapping of spin config- 
urations it is many-to-one: for a given choice of spin-1/2 
loops, there are many possible ways to occupy the re- 
maining edges of the lattice with spin or 1 states with- 
out violating the vertex condition. Second, the matrix 
elements of $p^i/2 are not all equal in magnitude, as is 
the case for the Z2 gauge theory. Third, matrix elements 
of the plaquette term in the Z2 gauge theory are also 
real and positive, while for SU(2)2 they can be negative 
or complex. 

It is thus somewhat surprising that the expectation val- 
ues of all operators required to describe the phase transi- 
tion are the same in both models. Essentially, this occurs 



because operators that commute with the conditions ( 53 ) 



can do only one of two things: either they measure the 
spin on each edge mod 1 (and are thus a linear combi- 
nation of transverse field operators), or they are linear 
combinations of $p,i (which in any case must act as the 
identity on our restricted Hilbert space) and i'p,i/2- As 
we show in Appendix [Cj the fact that the commutators 
of these operators have the same structure as their ana- 
logues in Z2 gauge theory ensures that their expectation 
values are identical in both models. 

For example, for F = any such expectation value can 
be expressed in terms the expectation value of an oper- 
ator that is diagonal in the the spin basis, after moving 
all $p_i/2, ^p.i terms to the right of all hg terms us- 
ing the commutation relations, and using the fact that 
^p,i/2|*o)r=o = ^p,i|*o)r=o = |*o)r=o- Expectation 
values of operators diagonal in the spin basis are insensi- 
tive to the relative phases between different spin config- 
urations in the ground states of the two models; they de- 
pend only on the probability of being in a given loop con- 
figuration. As we discuss in Appendix [Cl in the SU(2)2 
model adding (or removing) a closed loop of spin-1/2 



edges halves ( or doubles) the number of possible integer 
spin configurations compatible with the given choice of 
spin-1/2 loops. However, this is exactly compensated for 
by the fact that the coefficients in the action of the pla- 
quette projector multiply such configurations by a factor 
of \/2 (or l/\/2) relative to the state from which they 
were derived. Hence the probability to be in a given loop 
configuration, and hence the expectation value of any 
product of transverse field terms in the F = ground 
state, is the same in both models. 



1. Z2 Transitions m SU(2)k 

A; = 2 is a special case, in which the condensed phase 
is a Potts gauge theory, and the transition is also of the 
form described in Sect. |III A 1[ For larger values of k 
this is no longer the case. Here we briefly comment on 
the more general situation. 

First, we have already observed that we can always add 
a transverse field operator of the form hi'^''^-' = (— 1)^*". 
For transitions where this is the only transverse field, if 
we restrict our attention to states satisfying 



1 

V 



Lfc/2J 



P,s 



w^m (55) 



(where V is an appropriate normalization, such that the 
operator has eigenvalues and 1) the mapping ( |54| ) be- 
tween the low energy degrees of freedom in the SU(2)fe 
model and the Z2 theory remains valid. The rules for 
allowed combinations of angular momenta at a vertex 
ensure that edges of half-integer spin form closed loops, 
allowing us to map from closed 1/2-integer spin loops 
to closed integer spin loops in the Z2 gauge theory. Here 
again, because the only excitations relevant to the transi- 



tion are vortex loops created by h. 



(fc/2) 



which are Ising- 



like, as the operator squares to the identity - one can 



show that all correlation functions in the limit ( 55 ) are 



equal to their Z2 analogues. This is discussed in more 
detail in Appendix [C| 

The solvable Hamiltonian describing the physics of the 
condensed phase is simply 



1 - 



1 

V 



Lfc/2J 

E 

s = l 



P.s 



(56) 



The Hilbert space in the solvable limit contains only in- 
teger spins on the edges. For fc > 2, however, the vertex 
condition on these integer spins is not that of a Zp gauge 
theory: for example, if /c = 4, we have 



1x1=0+1+2 



(57) 



and vertices with 3 spin 1 edges, or two spin 1 edges and 
a spin 2 edge, are allowed. We will see in the next section 
how to determine the topological characteristics of these 
condensed phases. 
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VI. TOPOLOGICAL ORDER AND 
CONFINEMENT IN NON-ABELIAN 
WALKER- WANG MODELS 

Having investigated the confining phase transitions in 
Abelian and SU(2)2 Walker- Wang models, we now show 
that many of the results found follow naturally from the 
mathematical structures ('pre-modular categories') used 
to define Walker- Wang models. While the methods of 
this section will not furnish us with the details of the 
phase transition, they do tell us how topological order 
changes for any confining transition in a Walker- Wang 
model. Before setting out a general recipe, we summarise 
the basic properties of general WW models and their cor- 

for details). 



responding categories in Sec. VIA 



see 



Then in Sec. lyiBl we show that the results of Sec. HVBl 
and Sec. |V A| are simply described in terms of the M- 
matrix (closely related to the S'-matrix) of the categories 
corresponding to the CAWW and SU(2)fc Walker- Wang 
models respectively. Following this, we describe the con- 
densation in SU(2)fc for the cases k > 2 which did not 
yield to the methods of Sec. |V] 



A. General Walker- Wang models and the M-matrix 

General WW models have a Hilbert space consisting 
of a p state system on each edge of the lattice shown in 
Fig. [4] We denote the p possible states by the labels 
{ao, oi, . . . , ap_i}, where ag is the identity element; for 
the models of |III| these labels were the possible spins Sg = 
0, 1, ...p — 1 with Se = corresponding to ao- Each label 
in the set has a conjugate label also in the set, which 
we denote a.^. (For the abelian models, s = p — s, while 
for SU(2)fc, a = a). As we saw in Sec. Ill if a 7^ a, to 
define Bp and Qv we must include arrows on each edge 
to specify whether the operator measures a or a. The 
general Hamiltonian has the familiar form 



H = \j2a-Bp)-Mj2Qv-rj2k 



(m) 



(58) 



where hi™^ is defined below, and the operators obey the 
commutation relations in equation ( 19 1 and ( 20 ) . For 



now we tune the model to an exactly solvable point by 
setting r = 0. The Qv term ensures the ground state is 
a superposition of configurations for which only certain 
combinations of the p labels are allowed to meet at each 
vertex. The allowed combinations are those that satisfy 
the fusion rules of the category, as in Fig. |6] (a) - we call 
configurations that satisfy this condition at every vertex 
'string-net' configurations. 

The plaquette term Bp is subtler. Just as in the 
abelian case detailed in Eq. plaquette violations can 
be labeled by their 'fluxes', which are drawn from the 
same set {ao, ai, . . . , ap_i} as the edges (see Sec. HI A). 



These fluxes can be measured by a plaquette operator 
analogous to <tp, which raises or lowers the edge labels 



on the boundary of the plaquette P according to the 
rules of the category (see Refs. [15] and [22] for more de- 
tail). Since Bp and Qv commute, the result is that the 
ground state is a superposition of different string-net con- 
figurations, with coefficients that are related by the set 
of graphical rules in Fig. [6f^b)-(e). Only certain combi- 
nations of graphical rules in Fig. |6][a)-(e) are consistent 
with one another; when they are consistent, the resulting 
structure is called a premodular category. 

The p X p-sized 'Monodromy' or Af-matrix (which is 
essentially the modular S'-matrix with a different nor- 
malization) shown in Fig. [6] contains much of the useful 
information in the category. We will use this object to 
ascertain the topological order of phases on either side of 
the confining transition. Before doing so, we note that it 
satisfies'^ 



\Mab\ < 1 MaO = 1 Va e {ao,ai, 



}■ (59) 



In Ref. [551 we showed that the spectrum of point defects 
(and the associated topological order) can be deduced 
from the Ai"-matrix (or the closely related S'-matrix) of 
the category. Our main result was that point defects 
carrying charge j ^ are deconfined in the bulk if and 
only if Mij = 1 for all possible labels P^. As an example, 
we now use the AT-matrix to shed light on Walker- Wang 
models in the p = 2 models discussed in Sec. |III A"T] The 
p = 2 Potts gauge theory is a Walker- Wang model based 
on the so-called Zj^^ category, while the p — 2 CAWW 

model is based on a category Zj^^^-* . Both categories have 
edge labels in {0, 1}, and fusion rules specifying that an 
even number of label-1 edges must enter each vertex. The 
full properties of these categories is listed elsewher^^, 
but here we require just their Af-matrices 



Mt.c. = 




M- 



3DScm ^ 



1 1 

1 -1 



(60) 



In th e case of the p — 2 Potts gauge theory (or "3D toric 
code"'2SEll) ^ we see that both columns of M are filled 
with 1, so that particles carrying label 1 are deconfined. 
On the other hand the Af-matrix of the p = 2 CAWW 
model (also known as the "3DSem modeF'^^S]) has only 
its first column filled with 1, and so there are no non- 
trivial deconfined particles in the bulk of this model. In 
the less trivial non-abelian example studied in Sec. |V A[ 
the Walker- Wang model based on SU(2)2 has labels in 
{0, ^, 1} and an Af -matrix 



SU(2)2 




(61) 



Thus the SU {2)2 model from Sec. VA has no decon- 
fined bulk species since no column except the first is all 
unity. Having stated the correspondence between the M- 
matrix and the spectral properties of the corresponding 
Walker- Wang models, we now use it to understand the 
condensation transitions between different models. 
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FIG. 6. (a) Represents the vertices allowed by the category; 
the ground state of a Walker- Wang model will involve only 
these types of vertices. The diagrams in (b)-(d) serve two pur- 
poses. Firstly, they tell us the relative amplitudes of string- 
net configurations in the ground state e.g. row (c) tells us that 
configurations related by removing a closed loop carrying la- 
bel a occur with a relative factor of Aa in the ground state. 
Second, these diagrams provide a neat graphical mnemonic for 
the definitions of string operators. Note that there is a rule 
conjugate to (e) obtained by turning the over-crossing into an 
under-crossing on the left hand side, and sending R — >■ R* on 
the right hand side. 



B. Confining transitions from tlie M-matrix 



Ramping up F in Eq. ( 58 ) drives the system through 
a confining transition, where loops proliferate. In this 
section, we will describe the effect of this transition on 
the spectrum of the theory. To begin, pick an edge e, 
and suppose it is labelled j in the h direction. To create 
a small vortex loop of flux m encircling e in a right-hand 
sense with respect to n, define 



j) = M„ 



J) 



(62) 



In the simplest case, starting with the F = ground 
state, the state formed by vi™'' \ GS) has flux on every 
plaquette except those bordering e. If we set hi™^ = 
(we™'' + (ve^V)/^, then the resulting operator creates 
a superposition of flux m and flux —m loop around the 
edge e. Using the properties of the M-matrix in Eq. ( 59 ) 
it is clear that the theory in the F — > oo limit consists 



of a string- net with edge labels j satisfying Mmj = 1, 
although we are unable to say precisely how the system 
behaves for intermediate F. 

Starting with a category with labels {qq, oi, . . . , Up-i} 
and proliferating vortex loops of flux m, we can use the 
M-matrix to answer most questions about the topological 
order of the condensed phase. Firstly, the labels present 
in the condensed phase are precisely {ag, ai, . . . , a;_i} 
such that Maim = 1- The model obtained in the F — > oo 
limit is a new Walker- Wang model based on a cate- 
gory which has has the same graphical rules as the old 
category, except the rules are now restricted to labels 
{ao,ai, . . . ,ai-i} i.e., this is now a subcategory of the 
original category. In particular the new M-matrix is sim- 
ply the I X / matrix My = Ma-s.^ ■ Rather non-intuitively, 
if particles carrying label a are conflned for small F (i.e. 
Mam 7^ 1 for some a^) they can be deconfined in the con- 
densed (large F) phase, since it is possible that Ms^. = 1 
for the I labels at present in this phase. 



1. Abelian categories 



We now reproduce some of the results of Sec. |IVB| 
using the Af-matrix formalism. The CAWW models of 
interest have a p-state system on the edge of each lattice. 



and the transition in Sec. IV B was driven by proliferating 
vortex loops with an equal superposition of flux labels in 
{m, 2m, . . . ,p—m} (see Eq. (31 )). However, proliferating 
only the ±m vortex loops gives the same ground state in 
the F — ^ cx) limit (though the behavior at intermediate 
F/A may differ). To parallel more closely our treatment 
of the non-abelian models we will use this alternative 
form for the transverse field here, as it reproduces the 



topological orders found in Sec. IV B for the condensed 
phase. 

In the uncondensed phase, the M-matrix takes the 
form Mab = e±*2™''/P with a,b € {0,1, ... ,p - 1}, where 
we take ± for p even or odd respectively. If m labeled 
vortices are condensed, edges with label a develop a gap 
of ~ F(l — Re [Mam]) > 0, so the labels surviving on the 
edges of the lattice are precisely those satisfying Mam = 1 
i.e. {0,p/m, 2p/m, . . . , (to — l)p/m}. It follows readily 
that as F — ^ oo, the ground state of Eq. (58) is precisely 



that of Sec. IV in the same limit. Moreover, this ground 
state can be thought of as arising from a Walker- Wang 
model with the same graphical rules as in the uncon- 
densed phase, but restricted to the surviving labels. The 
effective Af-matrix of the new model is an to x to ma- 
trix of the form M^y = e±*27r:ryp/m^ ^ ^-^j^ rescaled labels 
x,y e {0, 1, . . . , to — 1}, and we can use this to diag- 
nose the topological order of the vortex phase. We now 
summarise the cases presented in Table |l] using this new 
formalism: 

• p/m^ = 2n: Clearly M^y = e±*4,.xyn ^ ^ 
x,y, so that all to remaining edge labels are de- 
confined. Therefore the model has the topological 
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order of a gauge theory. A pair of particles car- 
rying label X have mutual bosonic statistics because 
the i?-matrix (Fig. [6]Je)) of the category satisfies 

= gTi27ra;^« ^ ^ (g^g ApP.IeI. 



• p/m^ = 2n + 1: Again M^y = e±*27.xy(2n+i) ^ 
so the model appears to have the topological order 
of a gauge theory. However in this case a pair 
of particles with odd charges x exhibit fermionic 
statistics, a fact that follows from the form of the 
i?-matrix Rl% = e^F*'^^'" = (-1)^ (see App.|e|. 

• Tpjrr? — njm an irreducible fraction: In this case 
gcd(p/m,m) = 1, which implies that M^y = 
g±i27i-xyn/m ^-^^ ^^^^ equal to 1 when xy — 
mod TO, which in turn implies that all point par- 
ticles are confined in the bulk. Thus the resulting 
model is a CAWW , and only has surface (rather 
than bulk) topological order. 

• p/m^ n/b with 1 < 6 < m: Letting n/b be a frac- again find that M^k 
tion in lowest terms implies that M^y = e**^'^^^"/''. 
Clearly to/6 columns of M are filled with 1, which 
implies there are {m/b) — 1 (non-trivial) deconfined 
species, and to(1 — 1/6) confined species. Thus the 
condensed phase has both bulk topological order, 
and surface topological order. 



the phase has the 'surface topological order' of the 3D 
semion model i.e. the p = 2 CAWW model. 

Condense to = |, with k odd: In the case that k is 
odd, proliferating the | vortex loops leads to a string-net 
state with only the integer labels {0, 1, . . . , ^^^j- This 
is because Mjk < 1 unless j is an integer. Specifying 
integer particle labels a = x,b = y leads to an effective 
M-matrix: 



• / (2x+l 

1 



)(2y+i) 



k+2 



k+2 



Sin (i^) Sin (i^) 



(65) 



,w 



1}. This new Af -matrix has 



xy < I for all x,y > 0, and therefore all excitations 



where x,y € {0, 1 

are confined in the vortex phase. The condensed phase 
is described by a category called 80(3)^^211. 

Condense integer m — ^: In the case the k is even, we 
1 precisely when j £ {0,1, ... , |}, 

hence proliferating vortices with flux | leads to a string- 
net phase with precisely these integer labels. In this case, 
however, two columns of the new M-matrix (obtained by 



2. SU{2)k 



restricting Eq. ( 63 ) to integer labels) are formed entirely 
of I's. These columns correspond to the 0, |-labeled par- 
ticles, which are the only deconfined particles. We saw 
an example of this transition in the k = 2 case where 
we proliferated the label 1 vortex, and we were left with 
Z2°'**' topological order. 



Having seen that the results of Table |T] are reproduced 
precisely by the M-matrix formalism, we now investigate 
a non-abelian example. The category SU{2)}^ has fc -I- 1 
labels {0, 1, . . . , |}, and the Af-matrix takes the forrrpSl 




We now summarise the results of proliferating label to 
vortices. 

Condense half-integer to; For half-integer to, one can 
check that Almj < 1 unless j = 0. Therefore, prolifer- 
ating half-integer vortices m confines all nontrivial edge 
labels, leading to the trivial string-net with only the zero 
label in the F — >■ oo limit. 

Condense integer to with < to < |; In this case, 
Mjnj < 1 unless j = 0, |. Hence proliferating to label 
vortices gives a string- net phase with two labels 0, |. The 
effective M-matrix for these remaining particles takes the 
form: 

Therefore, in the case that k is even, the resulting phase 
has the topological order of a Z2 gauge theory (or equiv- 
alently the 3D toric code). In the case that k is odd. 



3. SU{N)k 

Among the possible condensation transitions in 
SU(2)fe, therefore, are two notable families: first, by con- 
densing an "integer spin" vortex loop, we arrive at a con- 
densed phase described by a Z2 model, which has bulk 
or surface topological order if k is even or odd respec- 
tively. Second, condensing the "highest spin" (fc/2) vor- 
tex loop generically results in a non-abelian topological 
phase which also has bulk (surface) topological order if 
k is even (odd). Moreover, this transition has the same 
low-energy description as the Z^"^*"^ transition in Sec. [h] 

Interestingly, an analogue of both of these exists in 
SU(N)fe models: proliferating vortex loops in the adjoint 
representation leaves a set of deconfined edge labels with 
Ztv fusion rules. Further, there is always an order N 
simple current; condensing the corresponding species of 
vortex loops produces a transition identical to that of 
the Z^°'*^ gauge theory. Though we will not derive these 
results here (see |42] for the necessary information about 
these categories), intuitively both families result from the 
fact that the group SU(A^) has a centre Z^r, which is also 
present in the related tensor category. The two families of 
transitions correspond either to condensing the vortices 
associated with this Zjy subgroup, or condensing vortices 
that have trivial Berry phase only with the Z^r subgroup. 
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VII. CONCLUSION 

In this work we have compared a family of phase tran- 
sitions in the relatively well-studied Zp Potts gauge the- 
ories with a related family in the recently introduced 
Walker- Wang models. In both models, the transitions 
that we consider can be understood as the condensation 
of vortex loops (i.e. loops of plaquette defects). Both 
admit an identical mathematical description of the tran- 
sition - allowing us to deduce from the work of Ref. |25] 
that all of the corresponding phase transitions are first 
order. For the Abelian models, "clock" -like variants of 
both Walker- Wang and lattice gauge models exist; in 
these models for p > 5 the single first-order transition 
splits into two second-order transitions separated by a 
gapless Coulomb phase^- . 

However, the relationship between the topological or- 
ders of phases connected by such condensation transi- 
tions is fundamentally different in the gauge theories and 
the Walker- Wang systems. The uncondensed (or decon- 

fPotts 



Potts 



z 

topological order ^ 



TP 



Trivial 



Zj°"^ topological order 



FIG. 7. (Color online) A sketch of the phase diagram for the 
fined) phase of the 2^°"'^ gauge theory is a topologically 2^"= model, as a function of the transverse field strengths Fi , 

{2tt vortex loops) and r2 {4n vortex loops), where charge is 
quantized in multiples of 1/4 (mod 1). For small Fi, there 
is a phase with Zj"**^ topological order. This is separated 
by first-order phase transitions from both the trivial phase at 
large Fi and a phase with Zg"**" topological order at large 
F2. For large Fi and F2, there is an additional first-order 
phase boundary between these two condensed phases. The 
three first-order lines meet at a triple point. The location of 
the phase boundaries is a rough guide, based on numerical 
analysis of a very similar model by Ref. 1431 



ordered phase with p deconfined charges (i.e. vertex ex- 
citations) - one for each possible value of the electric 
flux (i.e. edge spin label). Depending on the value of p, 
there may be several possible condensation transitions. 
For any p, we can simultaneously condense all magnetic 
fluxes. This confines all electric flux lines, thereby confin- 
ing all charges and completely destroying the topological 
order. If p is not prime, it is also possible to condense a 
subset of the magnetic fluxes. This confines only a subset 
of electric flux loops (corresponding to the charges that 
have nontrivial Berry phase with the condensed vortex 
lines), and leaves the remaining electric fluxes and their 
corresponding charges deconfined. Thus there are also 
condensation transitions between the Zp°"'' phase and a 
Z^°"'' phase, where m divides p. An example of the gen- 
eral structure of these phase diagrams is shown for p = 4 
in Fig. [7] In all of these transitions, there is a reduction 
in the ground state degeneracy, and topological entangle- 
ment entropy, as the system enters the condensed phase. 

In contrast, in the uncondensed phase of the CAWW 
models, there are no deconfined excitations in the bulk. 
The ground state is nonetheless a "string-net" , with loops 
(or nets, for p > 2) of non-trivial "electric flux" (i.e. non- 
zero spin, in the language of Sect.s [ll| and III). (Here 
we use "electric" and "magnetic" flux by analogy with 
the corresponding objects in the Zp°"'' gauge theory, al- 
though for the CAWW models the analogy is not ex- 
act). Because the ground state contains extended loops 
(or nets), here too we can add a perturbation that con- 
denses loops of "magnetic flux" (i.e. plaquette defects), 
driving a phase transition that confines "electric fiux" 
(i.e. edges carrying certain spin labels). As for the Z^"**** 
case, simultaneously condensing all magnetic flux loops 
engenders a transition to a trivial phase, in which all elec- 
tric flux lines are confined. If m divides p, however, we 
may once again condense a subset of the possible mag- 
netic fluxes, leaving m—1 deconfined types of electric flux 



loops. Surprisingly, many of the latter condensed phases 
have bulk topological order and deconfined bulk excita- 
tions: the topological ground state degeneracy grows as 
the system is driven into the condensed phase. 

As we discussed in Sect. |III[ in general we cannot rule 
out the possibility that the uncondensed phase of these 
CAWW models is in fact connected to the trivial phase 
in which all electric fluxes are confined, once we allow 
perturbations that introduce vertex violations into the 
ground state. (We note, however, that there are Abelian 
Walker- Wang models similar to those discussed here that 
are known to be symmetry-protected phases.^^ Interest- 
ingly, when a subset of the possible vortex loops is con- 
densed, there are transitions across which the topological 
order changes; these are necessarily transitions between 
distinct phases. Fig. [8] sketches the possible forms of the 
phase diagram for p = 4, as a function of the two trans- 
verse fields Fi,F2 which create vortex loops of magnetic 
fiux 27r and 47r, respectively. For small Fi,F2, there is a 
CAWW region. This region is separated from the regions 
at large Fi, F2 by a first-order phase transition. For large 
F2 , the system is in a Z2 phase (with the slight twist 
that the charges are fermions rather than bosons). Since 
this phase is topologically ordered, the phase boundary 
must persist in the presence of arbitrary perturbations to 
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FIG. 8. (Color online) A sketch of the phase diagram for the 
p = 4 CAWW model, as a function of the transverse field 
strengths Fi, (27r vortex loops) and r2 (47r vortex loops), 
where charge is quantized in multiples of 1/4 (mod 1). For 
small ri,r2, there are extended loops of electric flux, as in 
the solvable p = 4 CAWW model. This region is separated 
by first-order phase transitions from both the trivial region 
at large Fi and a phase described by a Z2 gauge theory with 
fermionic matter sources at large F2. As we have shown, the 
locations of the phase boundaries are identical to those of Fig. 



El 



the Hamiltonian. For large Fi, the system is in a trivial 
phase with all electric fluxes confined. The arguments 
of Appendix [C] ensure that the phase boundaries match 
those of the analogous Zp Potts-type model shown in Fig. 
[7j However, since neither topological order nor symme- 
try distinguish this phase from the uncondensed CAWW 
phase, the phase boundary need not persist when we al- 
low arbitrary perturbations (which do not commute with 
the vertex term Qv) to the Hamiltonian. 

We have also briefly discussed transitions in more com- 
plicated CAWW models, for which the surfaces admit 
states with non-Abelian topological order. In certain 
cases we are able to identify the nature of the phase tran- 
sition here too (which is again first order) by mapping the 
transitions onto transitions in the Z^"*'^ gauge theories. 
We have also outlined a completely general procedure for 
deducing the topological order of any condensed phase, 
and identified transitions from CAWW models with the 
surface topological order of SU{2)k Chern-Simons theo- 
ries to phases with bulk topological order. Thus the fact 
that a condensation transition that partially confines the 
allowed fiux (or spin) labels can deconfine some of the re- 
maining charges in the bulk is not unique to the Abelian 
case. 

Returning to our Zp°'*^ gauge theories, we note that 
in these systems there are generically two types of con- 
densation transitions: the first-order confining transition. 



in which vortex loops or monopoles proliferate, and the 
Higgs transition, in which charges condense. (The Higgs 
transition is second-order for p — 2, and first-order for 
p > iP^. It is known'^ that the Higgs and confined 
phases are not distinct (provided we condense the fun- 
damental charge or fundamental flux). The qualitative 
form of the phase diagram in this more general case is 
shown in Fig. |9j for the example p = 2. 

It is thus natural to ask about the analogue of the 
Higgs transition in the CAWW models. To drive such 
a transition, we must introduce an additional transverse 
fleld term (of the form r^We) which raises the value of 
the spin on each edge e. Unlike the case of 2^°**'' gauge 
theories, however, such an operator fails to commute not 
only with the vertex terms at each end of the edge, but 
also with the plaquette terms on some of the surround- 
ing plaquettes. Thus we can enter the Higgs phase only 
by making F™ large compared with both the coefficient 
of the vertex term and the coefficient of the plaquette 
term. Because of this complication, the Higgs transition 
in the CAWW models does not map onto that of the 
gPotts gauge theories in a straightforward way. Indeed, 
it is not obvious that there is a phase transition at all, 
given that there is no clear indicator that the CAWW 
model represents a distinct phase. What is clear, how- 
ever, is that for sufficiently large F™ the model is again 
in a trivial phase, with the spins on each edge diagonal 
in the We basis in the limit F^ — > 00. Further, as this 
trivial phase (together with the trivial confined phase) 
is identical to its counterpart in the 2^°"'^ gauge theory, 
the arguments of Ref. [T7| ensure that the F — > 00 and 
F,„ — 00 phases are connected. A qualitative sketch of 
the resulting phase diagram is shown in Fig. [9] 

Of course, in the Zp°"'' gauge theory, we may also con- 
dense charges other than the fundamental one, to obtain 
Higgs phases with non-trivial ( Z^°'*'') topological order. 
In the CAWW models, such transitions appear not to 
produce new topologically ordered phases, however. 

In summary, this work has detailed just some of the 
many possible phase transitions that can occur in the 
Walker- Wang models. Though in the limits that we are 
able to describe, the transitions themselves are not ex- 
otic, the phase diagram has some very surprising features, 
including the possibility of producing topological order 
from a phase which cannot be distinguished, by either 
topological order or symmetry, from the trivial phase. 
We hope that future numerical or analytical work will be 
able to shed light on this rather surprising phenomenon. 
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Appendix A: General form of the plaquette 
projector for Confined Abelian Walker- Wang models 

Here we give a g ene ral expression for the phase oper- 
ator Qp^rn in Eq. (21 1. As discussed in the main text. 
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FIG. 9. (Color online): A sketch of the phase diagram both Z2°^^" (a) and p — 2 CAWW (b) models, with an addional transvere 
field term Fm that can be used to drive a Higgs transition. For the Z2°^^^ model, there are two phases: the topologically ordered 
^Potts pjjg^gg^ g^jjfj ^jjg trivial phase which is obtained by condensing either charges or vortex loops. These are separa ted by a 
first-order transition (solid line) along the F axis, and a second-order (dotted line) 4D Ising transition along the Vm axis^^^. A 
first-order line emerges at the point where these two meet, terminating at a second-order critical point. The phase boundaries 
sketched here are based on the numerical work of Ref. 1451 For the CAWW model, we do not expect the phase diagram to 
match that of the Z2°**'^ model once F^ > 0, so that the loci and nature of the phase transitions are largely currently unknown. 
In particular, though the Higgs phase clearly exists (and has the same ground state, in the limit Tm — >■ oo, as for the Zf""** 
case), we do not know at present whether there is a phase transition along the Vm axis, or indeed to what values of Vm the 
transition along the F axis persists. 



0p,m is diagonal in the spin basis, and has eigenvalues 
e""™/P, < n < 2p — 1. Thus 6p,m can be expressed in 
terms of a product of the operators Se,± acting on certain 
edges proximate to the plaquette P. In order to express 
0p,m in this way, however, for p > 2 we must excercise 
some care about our choice of Se,-|- versus Se,_. To avoid 
confusion, our choice is indicated in Fig. 10 an arrow 



pointing in the x,y, or z direction on edge e indicates 
that we use Se_+ and Wg _|_; an arrow along —x, —y, or 
—z on e indicates that we use Se^_ and We,_. 

Let us begin by defining operators that measure the 
spin (rather than its exponential) on each edge: 



5e,+ pe/ — *ePe/ 



{p- Se)\Se) 



(Al) 



where G {0, 1, ...p — 1}. In terms of the operators S 
introduced in Eq. ([T]), we have 



Se,± = Exp [i (27r/p) Se,±] 



p odd 



(A2) 
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0+ 1 



FIG. 10. (Color online): Orientations of edges used in the 
definition of the Walker- Wang Plaquette operator, shown for 
the 3 different types of plaquettes on the point-split cubic lat- 
tice (see main text). Edges in dP are shown in bold; edges 
in *P are not bold. The definition of the plaquette operator 
requires 2 special edges in *P, shown here in blue and indi- 
cated with the letters O and U. For p even we also require 
two special edges O-l-l, U-l-1, which are shown here in red. 



(these are the non-bold edges in Fig. 10).) Each plaque- 



tte contains one "O" and one "U" edge, indicated in blue 
in the Figure. 

For p odd, we take: 



For technical reasons (which we discuss shorttly), the 
definition of the operator 0p^m depends on whether p 
is even or odd. To define 9p,m, we identify two special 
"crossed" edges in *P. (Recall that *P is the set of all 
edges connected to, but not bordering, the plaquette P 



Exp 



.TT to(p- 1) „ , 

* 1^0,1^0 — ^V,vu, 



(A3) 



where the subscripts U (O ) denote the U (O) edges, with 
orientations chosen as shown in the Figure. 
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p even 



For even p, the choice of Op,™ in Eq. (A3) does not 
lead to a well-defined Walker- Wang modeL This is be- 
cause changing the spin on an O or U edge by p leads to 
an overall phase of Exp [itt m{p~ 1)] = (—1)™. This is 
inconsistent with the notion that the spin state on each 
edge is defined only modulo p. 

Instead, the phase operator that we use for even p can 
be expressed as: 



6Rm = (-l)"Exp 



.TT m \ ^ 



P 



yp.m 



(A4) 



where the values of Vi = ±1 are indicated in Fig. [TOj as 
explained above, and flp.m is defined below. Note that 
X^e G*P is always a multiple of p, provided that the 
vertex condition is always satisfied. (When the vertex 
condition is not satisfied, the expression for Qp,m is more 
complicated, but can be deduced from the form of the 
Hamiltonian given in Ref. 1151 ) Thus for m even, the 
first two terms in the product ( |A4[ ) are always 1. For m 
odd, their combination is negative if the net spin of all 
edges leaving the plaquette (with the orientations show 
in Fig. 



an odd multiple of p. 

With the same definitions of O and U edges as used 
above, we have: 



?p,m = Exp 



TT m 



(A5) 



(sign(su+i,^u+i 
(sign(so+i,,.o+i 



ml)**"'"" X 



where U + 1 (0 + 1) are edges in dP that are adjacent 
to the O and U edges, as shown in the Figure. Note that 



here we take sign(O) = 1. 

Evidently, since p—1 and p are coprime, in either case 
the possible eigenvalues of Qp^m are e''^"™/?', Q < n < 2p. 

Readers familiar with TQFTs may wish to note that 
for p odd, we use the category SU{p)i (or zj^^^^^^^\ 



in the notation of Ref. [38|) . while for p even, we use 
U{l)p/2 (or zj^^^^).E3 For odd p the 6j symbols are all 
unity, and Qp^m is defined purely by the R matrices. 
For even p the 6j symbols may be n egative; the term 

Exp 



■Trm 



in Eq. (A4) gives the net ef- 
fect of the 6j symbols required to ruse the spin m with 
the spins on the edges of P, with orientations given in 
Fig. flOj (This term is also present in the 2D Levin- 
Wen^^ models, on which Walker and Wang's construc- 
tion is heavily based.) The term 9p_m is unique to 3D, 
and accounts for crossings of the string labelled m with 
certain edges in *P. (This is explained in more detail in 
Ref. HH). The choice of "O" and "U" edges depends on 
the projection (or angle of view) used to determine these 
crossings. 

3. Spectrum of the plaquette term 



10 1 is an even multiple of p, and positive if it is q 



Let us begin by deducing a few key properties of 
Qp^m necessary for understanding the 



P,m 



($p) 



plaquette projector. First, we prove that 
Bp^m — {j3p^ij 



P,rn 



(A6) 



{Qp,r. 



For odd p, this is immediate since 
acts only on edges that are not raised by ^p), and from 
Eq. (A3) it is clear that Qp^m = {^p,raj ■ For even p. 



it is a question of verifying that the two sign terms are 
indeed shifted appropriately. (All other terms commute) . 
To see this, observe that: 



(sign(su+i,,.u+i - m)) $p|su+i) $p sign (su+i,,.u+i - mod (to + l,p)) sign (su+i,,.u+i - l) |su+i) 



(A7) 



This holds because (with the choice of orientations given in Fig. 10) <i>p lowers the eigenvalue of su+i,i/u^.i by l(mod 
p)\ to compensate for this we have raised m by 1 (mod p) on the right-hand side. (Evidently, this identity also holds 
if we replace U by O). From this, we deduce that: 



$P (sign(su+i,i.„+i - (sign(so+i,,.o+i - 1))' 



(sign(s 



U+l,i/u+i 



I))"'"'"" (sign(so+i,^o+i -m)y-"° 



(A8) 



as required. 



Second, it follows that 



This is because {^p^ = 0p,p = 1- For p odd, the latter 
follows immediately from Eq. ( A3 ) , as p — 1 is even. For 
{Bp-if = ($p)^0p,p = 1 (A9) P <5ven, we have that Exp i7r]X,e.e*p ^e.,...] = 1' ^^^^ 
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(sign(so+i,.-o+i))^°'''° = 1- 



4. Confinement and deconfinement in Abelian 
Walker- Wang models 

We next show in more detail how the presence of 6 in 
the plaquette term leads to confinement of all excitations 
in the bulk and why eliminating certain labels from the 
effective action can produce deconfined excitations. 

There are two types of excitations in these models: 
vortex loops (which always incur an energy cost per unit 
length, for F = 0) and vertex defects, which correspond 
to deconfined charges in the Potts gauge theories. A 
pair of vertex defects at vertices vi and V2 is created by 
raising the spin along a series C12 of edges connecting 
vi to V2- We will call the operator that does this V^_^^. 
The excitations are deconfined if and only if there exists a 
raising operator that commutes with Bp for all P (except 
possibly those plaquettes of which vi and V2 are vertices). 
This is most easily seen at the exactly solvable point T = 
0, where the eigenvalue of Bp is conserved. Since the 
spectrum of Bp is discrete, such strings either cost no 
energy (if they commute with Bp) or incur a fixed energy 
penalty per unit length (if they do not) . 

Consider the effect of Vq_^^ for a trajectory C12 that 
contains, say, the O edge (but not the U edge) of some 
plaquette P. If Vq^^.^ simply raises the spins on all edges 
by r, then we have: 



-inm{p-l)/p ry^^^^^Qp^^ 



p odd 
p even 
(AlO) 

The commutators are unchanged if C12 contains the U 
edge (but not the O edge), since the O and U edges by 
definition have opposite orientations. 

We could also consider operators that simultaneously 
raise the spin on edges in C12, and multiply the result 
by a phase depending on the spins on edges in *Ci2 (i.e. 
on the set of edges connected to, but not on, the path 
Ci2). For example, if our operator acts on the O edge of 
a plaquette P, we can eliminate the phase difference be- 
tween BpmV( 



t 



and VQ^^Bp^m by adding to a phase 



of the form Exp 



for p odd (or 



-SO+l,vo+iSO,vo 

its equivalent for p even) . This will give an operator that 
commutes with Bp, but fails to commute with the pla- 
quette term on at least one of the neighboring plaquettes 
(with which it commuted previously) . Thus in general, it 
is not possible to create an operator Vci2 that does not 
incur a finite energy cost per unit length. 

Next, let us consider what happens when we restrict 
ourselves to a subset of the possible values of m. Suppose, 
for example, that both m and r must be multiples of p/s, 
where s\p -i.e. m = ap/s,r = bp/s. With this restric- 
tion, the possible phases appearing in the commutator 



of 0p,m and Vp^^.^ arc integer multiples of mhp/s^. If 
bp/s^ is not an integer, then the situation is as described 
above: if C12 contains both the O and U edges of P, 

-bp/s ' ®^'P/s 7^ ^' there is no way to adjust W 
that avoids paying a finite energy cost per unit length. 
Essentially the same reasoning can be used to show that 
there is no operators Wc that raises all spins on the non- 
contractible curve C hy r = bp/s and commute with the 
Hamiltonian. 

Conversely, if bp/s"^ is an integer, then there is an oper- 



ator Vc^^-r that commutes with Bp ap/s for every integer 
a and plaquette P. If, for C12 containing either the O 

or the U edge of P (but not both), V^^^ r' ^p,p/s — 0) 
then simply raises all spins along C12 by r, as for 



the Potts gauge theory. If jv^t^ 



P,p/s 



0, then 



Ci2;r 



also contains phases depending on the spin la- 
bels in *Ci2; nonetheless it is possible to assign these 

0. A simple way to see that 



such that 



there must be deconfined particles in this case is to note 
that y}j^^ commutes with lip, ap/s whenever C12 contains 
both the O and the U edge of P, allowing us to separate 
charges along certain directions without any energetic 
penalty. By adding extra phases to Vc^^i we can allow 
separations in any direction. 

To summarize, there is one deconfined particle species 
for each b such that e^'^^vh = as claimed in the 
main text. For each deconfined particle species, there is 
a corresponding operator Wc (for each non-contractible 
curve C) that commutes with the Hamiltonian, and can 
be used to map between topologically distinct ground- 
state sectors. 



Appendix B: Counting spin configurations in general 
Potts gauge transitions 

Here we verify that the condensation transition involv- 
ing he"*^ in a Zp°**** model can be mapped exactly onto 

a transition involving he^'' in a 2^°^** model. 

Our objective is to show that, for given F, A, the expec- 
tation value of any operator involving the flux excitations 
present in the ground state (i.e. those created by hi™'') 
is the same as the expectation of the analogous operator 
in the model. (Since the excitations that condense 

in both cases arc loops of flux l-Knmjp, the analogous 
operator is simply the operator creating the same config- 
uration of flux loops). 

Intuitively, this follows from the fact that each inser- 
tion of { S'J 



creates one flux loop (with flux 0p = 

2'!T{m/p)n) encircling the edge e, and that the rules for 
combining the flux loops are the same in both models. 
This ensures that the number of ways of obtaining a given 
conflguration of flux loops within perturbation theory is 



25 



the same in both cases. 

Let us see how this arises at the level of configura- 

tions. Since l^'p) commutes with all operators in 

the Hamiltonian, energetically there is no difference be- 
tween an edge of spin s and an edge of spin 'pjm + s 
(mod p). In other words, the many-to-one mapping be- 
tween configurations of the Zp""** model and the Z^""** 
model 



Se + np/m ( mod p) — ?> Se 



(Bl) 



preserves the energy of any spin configuration. In par- 
ticular, it does not affect the configuration of those vor- 
tex loops created by hi™'', as ^^™'' 



, (We,±) 



= 0. 



Hence the expectation value of any operator involving 
the vortex loops (i.e. any operator that is related to the 
degrees of freedom involved in the phase transition) is 
unaffected by this mapping. It follows that the phase 
transition must be identical in both systems. 

Because the mapping is many-to-one, it fails to cap- 
ture the physics of the degrees of freedom that are not 
involved in the phase transition. Indeed, the mapping 
sends m distinct spins of the Zp°"*' model to a single 
spin in the 2^°"'' model. Hence there are m deconfined 
spin variables in the condensed phase, as opposed to a 
single (sg — 0) deconfined spin variable in the condensed 
phase for the 1^^°^ case. 



Appendix C: Equivalence of Phase transitions in 
Potts and Walker- Wang models 

Here we give detailed arguments as to why certain 
phase transitions in Walker- Wang models must be of the 
same type as transitions known in the Zp°**^ models. 



1. Abelian Walker- Wang models 

We begin with the Abelian Walker- Wang models. Here 
all possible phase transitions can be mapped onto tran- 
sitions in the 2^°"*^ models. In practice, the arguments 
presented here do not depend on the relative coefficients 
of the different powers of Bp i and Se.+ in the definitions 
of Bp and hg, and can equally be applied to variants of 



the CAWW models, such as the Zj, lattice gauge theory 
of Sect. 



HC 



or the Fi — r2 models discussed in the 

Conclusion. 

Specifically, we will show that the expectation value of 
any operator composed of an arbitrary sum of products 
of flux measurements and spin measurements is identical 
in the two models, for any value of F/A. In the CAWW 
models, flux measurements are carried out by the oper- 
ator Bp,i = {^pQp^ij , whose eigenvalues we identified 

with the fiux through a given plaquette. In the Zp°**'' 
model, a flux measurement is given by Bp i = <^p. The 
two flux operators differ by a set of spin measurements, 
carried out by Se in both models. Hence the set of opera- 
tors that we allow is identical in both models; however, it 
is convenient to identify our "flux" operator as the object 
that commutes with the Hamiltonian for F = 0. 



First, we note that since 



= 0, we have 



^"/^Se,, f<lpep,i) eedP 



otherwise 
(CI) 



In the Zp°**'' models, the identical relationship holds, 
with ^$pOp,i^ replaced by ^p. In either model, there- 
fore, any product O of Bp.,i and Sej,,,, can be expressed 



in the form: 



a = e-2.«/p-Q|^S^ j 'l[(^Bp^ ,^ ' (C2) 



where q is the same integer for the Zp°"'' and CAWW 
cases. 

Second, observe that in both models, the exact ground 
state 1^*0) for any value of F/A can be expressed in the 
form: 

1*0) -E"""^^'^ n heJ*o)r=o (C3) 

a fee{c„} 

where {Ca} is an arbitrary collection of edges in the lat- 
tice, and a„ are arbitrary coefficients. This is true every- 
where in the phase diagram, though the coefficients Cq, 
can be evaluated perturbatively only in the uncondensed 
phase. Thus the expectation value of an operator O can 
always be evaluated via: 




If O is a sum of products of Bp^ and Se, we may use Eq. (CI) repeatedly to move all Bp^ operators to the right. 
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arriving at the form: 



(a)r.A = (*o|r=oEk'/'0 n he J I n h.J5:cAl*o)r=o 



(C5) 



a, 13 



vfee{c„} 



where Oj is of the form (C2), and the coefficients Cy are necessarily the same in both models as they follow directly 



from Eq. (CI). Since -Bp,i|^'o)r=o = |5'o)r=o, we have 



(^^>r,A = E (« 



(r/A)V„(r/A) 



r-0 



n n h,in(se 



o/r=o 



(C6) 



In other words, we can express (O) as a sum of expecta- 
tions of operators that are diagonal in the spin eigenba- 
sis, with coefficients that are completely determined by 
{a^a^^^} and the relation (CI). 

To complete the proof, we must show two things. First, 
that if Os is an operator diagonal in the spin basis, 
(^o|r=oC^s|^o)r=o is the same for both mo dels; and sec- 
ond, that the coefficients 02"^"^ in Eq. (C3) are the same 
for the CAWW and Z^""" ground states. 

For the first item, we observe that the F = ground 
state can be constructed from the trivial state |0) = 
rie = 0), via 



I* 



o/r=o 



n^pio) 



Since Bp — Bp (or more generally, BpiBp 



(C7) 



Bp, which 



follows from Eq. (A9)), this is an eigenstate of the pla- 



quette projector with eigenvalue 1; it also has no vertex 
violations, since Bp and Qv commute. Hence 



(^'o|r=oa|*o)r. 



(C8) 

Since by definition the spin is on every edge in the triv- 
ial state |0), it follows that this expectation value is com- 
pletely determined by the coefficients that result from 
moving Og past the product over plaquette projectors. 
These coefficients are necessarily the same in both mod- 
els (see Eq. (CI)). Hence this expectation value must 



be identical in the CAWW and Z^"**'' ground states. In- 
tuitively, this is because the matrix elements of $p and 
(^^pOp^ij between any pair of spin configurations differ 

only by a phase; since different spin configurations are 
orthogonal, operators diagonal in the spin basis cannot 
measure these phases. 

Notice that in the Z^"***^ case, Eq. (C7| describes 



only one of the possible ground state sectors in periodic 
boundary conditions. Hence any statements that we can 
make about the correspondence in the phase transition 
between the two models are valid in the trivial ground 



state only. This will not affect the thermodynamic prop- 
erties of the transition, however. 

Finally, we must show that the coefficients a^a^^'' are 
the same for both models. This is clearly true to any 
order in perturbation theory: at F = the spectra of the 
two Hamiltonians (in the topologically trivial sector) are 
identical, and the perturbing term F hg has the same 
effect on both of these spectra. Beyond the phase transi- 
tion perturbation theory is no longer valid; however, we 
may re-express the wave-function (|C3|) as: 



" A;e{C„} 



(C9) 



where {Ca} is a set of plaquettes. The coefficients ^L^^^"* 
must be equal in both models to any order in perturba- 
tion theory, for the same reasons as in the uncondensed 
phase. Further, if e is an edge of P, in the Z^"**** case we 
have: 

p— 1 p—i 

^ n=0 m=0 



p. 



Bp^i 



(CIO) 



(If e is not an edge of P, the two operators commute). 
Hence 



p— 1 p— 1 



pi27rmn/p|r|\ _ 



1, 



|0) = -|0) 

p 

(Cll) 

Thus the wave function ( |C9| can also be expressed as: 

E^i'^"^ n BpjO)=J2b^^/^^ n {pK)\^o)r=o 
" feefCa} " fcefc*} 

(C12) 

where {C* } contains at least one edge of every plaquette 
not in {Ca} (and no edges of the plaquettes in {Cq}). 
(Including multiple edges of a particular plaquette will 

not affect the outcome). Since bi^^^^ are equal in both 
models, this proves the desired result. 
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If he OT Bp are not of the Potts form (in the sense that 
they do not contain all powers of S or $ with equal am- 



plitudes) Ens, (CIO) and (Clll no longer hold. However, 
in Eq. (C3| we need not restrict he to the transverse 



field operator in the Hamiltonian in order for the result 
to be valid; it suffices to express the ground-state wave 
function in terms of operators diagonal in the spin basis 
acting on the F = ground state. We therefore replace 
Eq. (|C12|) with 



J2bi'^^^ n n (phe,)|vI/o)r=o (C13) 

Here C* is defined as above, and Ca contains the re- 
maining edges of the lattice (notably, it contains at 
least one edge of every plaquette in Ca). Here he — 

Sn=o (^'^'+) ' ^^^^ chosen so that that 



P 



heE^p"l|0)=^P|0) 



(C14) 



This gives an expression of the required form for the 
ground state at large F/A, with coefhcients that are nec- 
essarily the same in both models. 



2. SU(2)fe Walker- Wang models 

Next we consider the set of transitions for SU(2)fc in 
which spin fc/2 vortex loops condense. The arguments 
we used in the Abelian case can essentially be applied 
here, as they depend only on the commutation relations 
between Bp i and the transverse field term, with a few 
modifications which we describe below. 

Specifically, if we shift the transverse field term by a 
constant, taking 



(-1) 



2Se 



(C15) 



then Se and $p.i/2 anti-commute, satisfying Eq. (CI) 
with p = 2. Hence Eq.'s (C3) - (C6) apply, and we need 



(r/A) 



only show that the expectation values of any product of 
transverse field terms in the unpe rturbed ground state, 
and the coefhcients a, 
of the Z2 gauge theory. 



in ( C3 1 arc identical to those 



For SU(2)fe, Eq. (C7) must be replaced by: 



I* 



S(7(2)fe> 



;r=o — 



f + $P 1 - V2$ 



Pa/2 



|0) (C16) 



The transverse field term (—1)^^= commutes with 1 -t- 
<i?p_i, but anti-commutes with $p.i/2. Thus moving the 



transverse field past the plaquette operators has the iden- 
tical effect as in the Z2 gauge theory, where 



I* 



/r=o 



p 



(C17) 



and Se commutes with 1, and anti-commutes with $p. 
It follows, by the same reasoning as used in the Abelian 
case, that 



(C18) 



where Q^^^^^^^'' measures spin mod 1 in the SU(2)fc 

model, and Ol ^' measures the Ising spin of the Z2 ga uge 
theory. Similarly the arguments leading to Eq. ( C12[ ) re- 
main valid, provided that we replace the trivial state |0) 
with the state obeying (53). This completes the proof. 



We emphasize that Eq. (C18) does not hold if we 



include operators that can differentiate between integer 
spins in the SU(2)fe models. Such operators will not com- 
mute with <I>p^i, and this component of the plaquette 
term (which is analogous to the identity component of 
the plaquette term in the Z2 case) is not left invariant by 
moving the spin- measuring operators to the right. 

As mentioned in the main text, this result is somewhat 
surprising in light of the fact that the matrix elements of 
^p,i/2 in the SU(2)fc model and $p in the Z2 model differ 
both in sign and magnitude, and because a given config- 
uration of half-integer spin loops in the SU(2)fc model 
corresponds to many different spin configurations. It is 
instructive to understand why these differences do not af- 
fect the expectation values of any operator that satisfies 
the conditions (53) in the F = ground states. This is 



most apparent in the SU(2)2 case, though the counting 
carries through for any k.^^ 

As for the Abelian case described above, operators sat- 
isfying ( 53 1 are insensitive to the phases of the matrix 
elements. This is clearly true for any operator diago- 
nal in the spin basis, and also (by the arguments above) 
for <I'p.i/2. Less trivially, the difference in magnitude 
and multiplicity of configurations turn out to cancel each 
other: in a given spin configuration in the SU(2)2 model, 
any contractible spin-1/2 loop can be removed and re- 
placed with one of exactly two configurations of spins 
and 1. This is because (1) the number of spin I edges 
external to any spin 1/2 loop must be even if Eq. (53) is 
satisfied; and (2) for a given loop on which the total num- 
ber of external spin 1 variables is even, there are exactly 
two configurations of integer spins on this loop which sat- 
isfy the vertex condition everywhere. The first condition 
follows from the fact (explained in more detail in Ref. I22|) 
that for F = the ground state of a Walker- Wang model 
contains only configurations that are "allowed" as link 
diagrams within the category; a diagram with a spin-1/2 
loop and an odd number of spin-1 external edges is not 
such a diagram. To see the second fact, we first note 
that if all edges external to the loop have spin 0, then we 
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may replace the spin- 1/2 loop with either a closed loop 
of spin or of spin 1. If there are at least two external 
edges of spin 1 , then the two ways to assign integer spins 
to the loop correspond to beginning at one of the spin-1 
external legs, and assigning a spin of 1 to one of the two 
edges on the loop with which it shares a vertex. After 
this initial choice there is no more freedom in the spin as- 
signments, provided that the spins of all external edges 
are fixed and we do not introduce any vertex violations. 

By the above reasoning, if we add (or remove) a closed 
loop if spin-1/2 edges, we halve ( or double) the number 
of possible integer spin configurations compatible with 
the given choice of spin-1/2 loops. However, this is ex- 
actly compensated for by the fact that the coefficients in 
the action of the plaquette projector multiply such con- 
figurations by a factor of \/2 (or l/-\/2) relative to the 
state from which they were derived. Hence the probabil- 
ity to be in a given loop configuration is the same in both 
models. Hence from the point of view of operators only 
sensitive to the spin mod 1, the ground-state loop config- 
urations of the SU(2)2 Walker- Wang model look exactly 
like the loop configurations in the ground state of the Z2 
model. 

Needless to say, here again the arguments apply only 
to the trivial ground-state sector of the Z2 theory; the 
other ground-state sectors (where they exist) do not have 
analogues in the Walker- Wang case. 



Condensing vortices with label m out of the Walker- 
Wang model leads to a phase described by a Zf^^ 
Walker- Wang model, where g — gcd{2qm,p). To see 
why, notice that the labels remaining in the condensed 
phase are precisely those a satisfying e'^'^'Ja™/p — i xhe 
most general values of a satisfying this equation comprise 
the set {0,1,2, g — 1} X p/g, and so precisely g labels 
remain in the condensed phase. 

The new model inherits the old fusion and braiding 
rules, as well as the M-matrix from Z^, but restricted 
to the g new labels. Indeed, the new model can be 
thought of as being based on the category Zf^^, which 
has labels x G {0,1, . . . , g — 1} related to the old la- 
bels by a = xp/g. We can understand the topological 
properties of this phase by examining the new M-matrix 
which takes the form M^y = e'^^'^^^/f. While the un- 
condensed model had c — gcd{2q,p) deconfined labels, 
the condensed model has c = gcd{2qp/g,g) of them. In 
the Zp gauge theories (which have q = 0), c = p while 
c = g. In the CAWW models, there is one deconfined 
label (a — 0), while c = gcd{p/g,g) for p even, and 
gcd{p{p- l)/g,g) = gcd{p/g,g) for p odd. 



Appendix E: Fermions versus bosons in 
Walker- Wang models 



Appendix D: More general Abelian Walker- Wang 
models and their transitions 



In Sec. |III| we considered CAWW models, which com- 
prise only a subset of all abelian Walker- Wang models. In 
this section we examine vortex transitions in more gen- 
eral abelian Walker- Wang models. Since the arguments 
of Appendix [C] can be used to map the phase transitions 
onto those of discrete Abelian gauge theories, we will fo- 
cus here on the nature of the condensed phase. It will 
be useful to speak about Walker- Wang models in terms 

of their corresponding categories. The p-state CAWW 

1/2 

should be thought of as being based on the category Zp 
for p even, and ^^^^ for p odd^® . Other examples in- 
clude the abelian Potts gauge theories of Sec. |ITj which 
can be viewed as Walker- Wang models based on Z^. 

In general, p-state abelian Walker- Wang models are 
based on categories Z^, where 0<q<p— lis either 
integer or half-integer, and where q can only ever be half 
integer if p is even. These categories are described in Ref. 
|38l The labels in Z? are {0, 1, ... ,p - 1}. The general 
M-matrix of category Z^ takes the form Mab = e'''^'"^°-''/P , 
and as before we can use this to understand the topo- 
logical order of the corresponding Walker- Wang model - 
remember that if M has c columns filled with 1, then 
the corresponding Walker- Wang model has c deconfined 
species, and a topological degeneracy of on the 3-torus. 
c = gcd{2q,p) in the case of Z^. 



In Sec. IV B we claimed that spin labels a for which 
the eigenvalue of 0p_a can be only ±1 are associated with 
fermionic "charge" excitations, while if = 1 then the 
associated charge is a boson. In the language of cat- 
egories this is the statement that Walker- Wang models 
based on Zp^^ (which can only be defined for even p) have 
deconfined fermionic and bosonic excitations, while those 
based on Z^ have all their point-like excitations decon- 
fined and bosonic. In this section we show how fermionic 
statistics arise in Walker- Wang models based on Zp^^. 

To begin, act on the ground state with a deconfined 
string operator which connects vertices 1 to 1, 

and creates two defects with conjugate charges a, a e 
{0, 1, . . . ,p — 1} (note a ^ a in general). Similarly use 
Wq^ _ I GS) to create an identical pair of defects of charge 
a, a at 2 and 2 respectively, thus producing the excited 



state W^^_W^^_ I GS) (see Fig. 11 a)). The prescrip- 
tion for forming string operators is discussed in detail in 
Ref In the cases of Zp''^ and Z^, the string operator 
creates deconfined excitations because it commutes with 
the plaquette operators along its length; this is possible 



cd 



1. 



because Zp''^ has a trivial M-matri:J2ll ]^ 

The exchange operation is only meaningful for iden- 
tical particles, thus we begin by exchanging sources 1 
and 2 (which have equal charges) which corresponds with 
adiabatically evolving the state Fig. 11 'a) to the one in 
Fig. [TT|^b) . This new state (a) will be equivalent to (b) , 
except for an overall phase. To evaluate this phase, we 
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state. 



FIG. 11. (a) shows point sources in an abelian Walker-Wang 
model with "Fermionic excitations", where 1, 2 have charge a 
and 1, 2 have charge —a. The strings represent the operators 
that act on the ground state to form these point excitations, 
(b) Shows the diagram resulting from the exchange of the 
defects. Plaquettes operators can be used to hop the string 
to the configuration in (c). The over-crossing in (c) can be 
removed at the cost of introducing a possible minus sign, and 
thus the diagram can be deformed back to (a) but with a 
phase of (—1)", which gives a fermionic or bosonic statistic 
contingent on whether a is odd or even respectively. 



Notice that the state in Fig.[TTJ^c) has an over-crossing 
of the string operators (from the point of view of the 
projection). The overcrossing can be remo ved and we 

d), but 
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can connect 1 to 2 and 2 to 1 to form Fig. 
there will be a phase associated with this process, as 
illustrated by 



= (-ir 



r 




„ „ (E2) 
which shows that ($p)"6p.a introduces a phase if it 
changes the number of over-crossings in a configuration. 
Notice that ($p)°6p.a = 1 on the ground state, so con- 
figurations that are related by removing such an over 
crossing have a relative sign of (—1)° in the ground state 
superposition. One can extend this reasoning to show 
that the state in Fig. 11 c) is the same as Fig. [lljd), 
except for this relative phase. Thus, the exchange of the 
two defects has resulted in a configuration identical to 
that in Fig. 11 |a), but with a phase of (—1)". 

One can more carefully derive this exchange phase 
by noting that, when a category has trivial M-matrix 
like Zp^^, the strings representing operators can be de- 
formed according to the graphical rules of the category 
seen in Fig. [61 This fact follows from the 'handle-slide' 



propert y ' of Walker- Wang models. To get between 
Fig. [njc) and (d) , one may use the manipulations of the 
string operators 



first note that if P is a plaquette touching a path C, then 



(l>p)'^ep,„W^c" I GS) - WS, I GS) 



(El) 



where C is the same path as C except it takes a detour 
around plaquette P (this step relies on the fact that Zp^^ 
has trivial M matrix). We can now use the fact that 
Qp^a — 1, and the fact this operator commutes with W^, 
to shown that W^, and have the same effect on the 
ground state. Using this idea repeatedly, it is easy to 
show that Fig. [Il|(b) and Fig. [TT|(c) represent the same 




(-If 




.{-If 



.(E3) 

which correspond with a use of the rules in Fig. |6j 
with Fig. ^e) needed specifically for the second equal- 
ity, where we used the fact that = (-1)" for Z^^^ 
The statistic arises precisely because the strings connect- 
ing the point defect keep track of the exchange. 
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